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INTRODUCTION 


As  equipped,  the  power  system  simulator  (PSS)  is  designed  for 
simulating  electrical  transient  behavior  in  a power  system  rated  up  to 
0.5  MVA  with  up  to  20  transformers.  The  basic  technique  for  this  simula- 
tion is  to  eaiploy  analog  computation  methods  to  pre-program  each  power 
element  (transformer,  inductor,  resistor,  etc.)  into  a simulation  module. 
The  input  and  output  terminals  of  these  modules  will  be  interconnected 
according  to  a signal  flow  diagram  of  the  selected  power  system  to  form 
a simulated  power  system.  If  a simulated  transient  is  injected  into  an 
input  terminal  of  the  simulated  power  system,  different  transient  re- 
sponses will  appear  on  the  terminals  of  each  module  simultaneously. 

These  responses,  therefore,  can  be  recorded  for  use  in  transient  analy- 
sis. Besides  the  transient  simulation,  PSS  also  can  be  used  for  studying 
electrical  steady-state  problems  at  various  AC  frequencies.  The  simula- 
tor can  be  equipped  to  double  its  present  capability. 


GENERAL  DESCRIPTION 


As  shown  in  Figure  1,  the  PSS  consists  of  three  bays.  On  the  top 
of  the  middle  bay  are  three  cages,  each  of  which  can  contain  a maximum 
of  eight  modules.  A control  panel  is  located  directly  below.  The 
buttons  on  the  panel  control  the  operational  modes  of  the  PSS  such  as 
the  compute,  hold,  set,  and  reset  modes.  The  output  of  each  module  can 
be  read  through  an  output  reader,  located  at  the  center  of  the  middle 
bay,  by  pressing  down  the  proper  button  in  each  row  of  the  button  matrix 
(lower  left  corner  of  the  control  panel).  The  lower  part  of  the  middle 
bay  contains  voltage  reference  modules  that  provide  a stable  reference 


voltage  for  analog  computation  and  buffer  amplifier  modules  that  provide 
a buffered  output  to  recording  or  controlling  devices.  The  bays  on  the 
left  and  right  are  identical;  potentiometer  panels,  used  to  scale  dovm 
input  signals  to  the  modules  by  a designed  factor,  are  in  the  middle  of 
each  of  these  bays.  Besides  the  potentiometer  panel,  each  of  these  two 
bays  contain  four  cages;  each  cage  can  also  accommodate  eight  modules. 

In  general,  the  PSS  consists  of  three  kinds  of  programmable  modules: 
(1)  simulating,  load,  and  transformer  modules  for  simulating  the  power 
elements  in  a power  system;  (2)  general  analog  modules  for  generating 
transient  signals  and  simulating  dynamic  loads;  and  (3)  a function 
generator  module  for  simulating  nonlinear  analog  signals.  In  addition 
to  these  basic  modules,  PSS  also  has  several  voltage  reference  modules 
for  providing  a standard  reference  voltage  and  two  buffer  amplifier 
modules  for  amplifying  output  power  to  drive  recording  equipment. 

Totally,  the  PSS  has: 


Modules 

30  load 

24  transformer 
5 summer/ Integra tor 
4 multiplier 
2 comparator 
1 function  generator 


Modules 

1 10-V  reference 
1 1-V  reference 
1 100-mV  reference 
1 10-mV  reference 

1 1-mV  reference 

2 buffer  amplifier 


The  PSS  was  designed  to  operate  at  a center  frequency  of  60  Hz  with 
a full-scale  amplitude  of  ±10-V.  The  nominal  frequency  range  of  simula- 
tion is  from  6 to  600  Hz,  which  can  be  exceeded  but  with  some  loss  in 
accuracy. 
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PROGRAMMABLE  MODULE 

i 


Programmable  modules  play  a very  important  role  in  PSS.  The  input 
and  output  signals  to  these  modules  are  designed  in  volts  with  their 

magnitude  in  the  range  of  ±10  V.  However  these  signals  need  not  neces-  j 

i 

sarily  represent  the  input  or  output  voltage  of  a simulated  power  ele- 
ment, they  can  represent  any  physical  quantities  such  as  current,  power, 
speed,  torque,  etc.  Because  the  magnitude  of  these  input  and  output 
signals  is  limited  to  ±10  V,  the  magnitude  of  these  signals  is  designed 
to  be  proportional  to  the  real  physical  quantities. 

Simulating  modules  are  designed  to  simulate  the  most  common  power 
elements  in  the  electrical  power  system  such  as  transformers,  resistors, 
inductors  and  capacitors;  thus,  they  are  the  most  useful  modules  among 
the  programmable  modules  in  the  PSS.  Since  each  power  element  can  be 
represented  by  a standard  mathematic  module  (the  characteristic  transfer 
function  CTF) , each  simulating  module  is  based  on  a mathematic  module 
and  has  the  same  analog  CTF  built  into  its  integrated  circuit.  The 
performance  parameters  of  the  CTF  in  each  module  are  adjustable;  there- 
fore, if  a power  element's  electrical  CTF  is  known,  the  CTF  of  the 
simulating  module  can  be  adjusted  to  match  the  electrical  CTF  of  the 
simulated  power  element.  Once  these  two  CTFs  are  matched,  the  power 
element  is  simulated  analogically  by  the  module. 

Next  in  importance  to  the  simulating  modules  are  the  general  analog 
modules,  which  actually  are  conventional  analog  modules.  They  are  more 
flexible  in  use  and  easier  to  program  but  require  more  modules  to  com- 
pletely simulate  a power  element  than  are  necessary  with  the  simulating 
modules.  Thus,  the  general  analog  modules  are  used  for  simulating  a 
transient  source  or  a special  power  element  that  has  not  been  pre- 
programmed into  the  simulating  modules.  Each  programmable  module  is  de- 
scribed in  the  following  sections  of  this  report. 
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Load  Module 


This  module  consists  of  two  identical  circuits.  Each  circuit  has 
its  own  independent  input  and  output  jacks,  adjustable  performance 
parameters,  and  overload  indicator.  Therefore,  the  module  can  simulate 
two  individual  loads  simultaneously.  The  analog  components  of  each 
circuit  consists  of  a summer,  an  integrated  operational  amplifier,  and 
an  inverter.  The  CTF  of  this  module  is  built  into  the  integrated  circuit 
amplifier.  Figure  2 shows  a simplified  circuit  diagram  of  this  ampli- 
fier circuit.  The  CTF  of  this  module,  expressed  in  Laplace  transforma- 
tion form,  is 


e.  ■ B S + 1 


(1) 


where  e 

0 

e . 

1 

S 

A 

B 


output  signal  in  volts 

input  signal  in  volts 

Laplace  operator 

adjustable  performance  parameter 

adjustable  performance  parameter 


The  CTF,  expressed  in  the  form  of  Equation  1,  is  the  most  common 
expression  for  load  element  in  the  power  system.  Figure  3 shows  three 
kinds  of  load  having  a CTF  similar  to  that  of  Equation  1.  These  load 
elements  are  (1)  the  resistor  (Figure  3a),  (2)  a capacitor  parallel  with 
a resistor  (Figure  3b) , and  (3)  an  inductor  series  with  a resistor 
(Figure  3c) . The  electrical  CTF  of  the  resistor  can  be  expressed  in  two 
ways , 


i 


R 


and 


(2) 


(3) 
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where  V = voltage  across  the  simulated  resistor 
i = current  through  the  simulated  resistor 
R = resistance  of  the  simulated  elements 

For  matching  the  CTF  of  the  load  module  to  the  CTF  of  Equation  2,  it  is 
simple  to  define 

e^  = V,  voltage  across  the  simulated  resistor 
e^  = i,  current  through  the  simulated  resistor 

and  adjust  the  performance  parameters  to 

A = -R;  B = 0 

For  matching  the  CTF  of  the  load  module  to  that  of  Equation  3,  the 
variables  in  Equation  1 will  be  defined  as 

e^  = i,  current  through  the  simulated  resistor 
e^  = V,  voltage  across  the  simulated  resistor 

and  adjust  A and  B to 

A = B = 0 

The  electrical  CTF  of  a capacitor  in  parallel  with  a resistor  is 

V = 5 (4 

i R C S + 1 

where  C = capacitance  of  the  simulated  element 
R = resistance  of  the  simulated  element 
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so  that  the  matching  of  these  two  CTFs  can  be  obtained  by  defining  the 
variables  in  Equation  1 as 

= V,  voltage  across  the  simulated  element 
e^  = i,  current  through  the  simulated  element 

and  adjust  the  performance  parameters 


A = - A;  B = RC 

The  load  shown  in  Figure  3b  has  its  electrical  CTF: 

i = 1/R 

V (L/R)S  + 1 

where  L = inductance  of  the  simulated  element 

R = resistance  of  the  simulated  element 

In  this  case,  the  matching  of  CTFs  is  to  define 

= i,  current  through  the  simulated  element 
= V,  voltage  across  the  simulated  element 

The  performance  parameters  A and  B will  be  adjusted  to 


Transformer  Module 


(5) 


The  module  is  used  to  simulate  a single-phase  transformer.  The 
basic  structure  of  the  module  is  the  same  as  the  load  module,  but  the 
difference  is  that  the  CTF  of  the  transformer  module  is  formed  by  com- 
bining two  independent  circuits  of  the  load  module.  The  module  thus 
contains  four  summers,  two  integrated  (circuit)  operational  amplifiers. 
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[ 
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and  two  overload  indicators.  A typical  equivalent  circuit  of  a single- 
phase transformer,  shown  in  Figure  4,  is  a two-port  T-network  with  its 
four  terminal  variables:  Vj,  V2,  i^,  and  12-  Tht  input  variables  of 
the  module  were  chosen  as  and  V2;  and  the  output  variables  as  i^  and 
i2.  The  loop  equations  for  this  equivalent  circuit  are 


or 


^1  = (i  ® ^ 2^1  h ^ K s(i^  + i^) 

^2  = (TS"2^)i2  " ^S(i,.i2) 

'^l  ■ ^^2  ■ (t  ® Fg)  ^^1  ■ ^2^ 


V + V = 

Vi  T V2 


2^]  (h*  V 


(6) 

(7) 


Let 


L,.  = L + f 

t m 4 


Equations  6 and  7 become 

2 G 


^1  ' ^2  _ 

Vj  - V2  G L S + 1 


2 G 


h * "2  _ 

+ V2  " 4 G S + 1 


(8) 


(9) 


In  comparing  these  Equations  8 and  9 with  Equation  1,  one  finds  that  the 
equations  are  in  the  same  form,  therefore,  the  transformer  module's 
structure  is  similar  to  that  of  the  load  module.  The  signal  flow  diagram 
of  the  transformer  module  is  shown  in  Figure  5.  The  signals  would 
sometimes  have  to  be  scaled  down,  scaled  up,  or  inverted  to  avoid  over- 
load, to  reduce  the  computing  errors,  and  to  take  into  consideration  the 
operation  of  the  amplifier  circuits.  The  techniques  for  matching  the 
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CTF  of  the  transformer  module  to  the  electrical  CTF  of  a simulated 
transformer  are  the  same  as  those  described  in  the  preceding  load  module 
section;  i.e.,  adjusting  the  performance  parameters  to  match  the  CTFs. 

Summer/ Integrator  Module 

This  module  contains  two  integrators  and  two  summers,  each  one 
containing  its  own  independent  circuit  and,  thus,  being  capable  of 
operating  independently.  The  feedback  element  in  the  integrator  can  be 
chosen  as  either  one  of  three  built-in  elements;  i.e.,  one  of  two  differ- 
ent capacitors  or  a resistor  with  fixed  resistance.  If  the  feedback 
element  is  a capacitor,  the  device  is  an  integrator;  but  if  the  feedback 
element  is  a resistor,  the  device  is  a summer.  Besides  the  regular 
input  and  output  terminals,  the  integrator  has  an  additional  IC  terminal 
for  the  purpose  of  setting  initial  conditions  in  an  integrating  opera- 
tion. The  summer  in  this  module  sums  up  all  incoming  signals,  inverts 
the  sign,  and  transmits  the  resultant  signal  to  the  next  stage.  Thus, 
the  summer  can  be  used  as  an  inverter.  The  summer  has  five  input  termi- 
nals, three  1-gain  terminals  and  two  10-gain  terminals,  and  four  common 
output  terminals.  The  operational  mode  of  this  module  is  controlled  by 
mode  control  buttons  located  at  the  lower  center  of  the  control  panel. 

Multiplier  Module 

This  module  is  a dual  module  which  contains  two  independent  multi- 
plying circuits.  Each  of  the  circuits  consists  mainly  of  an  integrated 
circuit  block  with  three  input  terminals  and  one  output  terminal.  Using 
different  patching  techniques  to  these  terminals,  the  integrated  circuit 
block  can  function  as  a multiplier,  divider,  or  square-root  device. 

Figure  6a  shows  that  the  integrated  circuit  becomes  a multiplier  by 
connecting  terminals  e^  and  e^  together.  In  this  configuration  e^  and 
e^  are  independent  input  terminals.  The  resulting  output  signal  from 
terminal  e^  is  equal  to  the  product  of  the  two  signals  from  the  input 
terminals  divided  by  10. 


8 


The  connections  shovm  in  Figure  6b  converts  the  integrated  circuit 
block  to  a divider.  The  input  terminals  are  e and  e , and  the  output 
terminal  is  formed  by  connecting  the  terminals  e^  and  e^.  The  output 
signal  is  equal  to  the  signal  from  terminal  e^  divided  by  the  signal 
from  e^  and  multiplied  by  10. 

Figure  6c  shows  the  connection  diagram  that  converts  the  integrated 
circuit  into  a square-root  device.  In  this  configuration,  the  terminals 
e^,  e^  and  e^  are  combined  together  as  the  output  terminal;  the  input 
terminal  is  e^.  The  resulting  signal  from  terminal  e^  is  equal  to  the 
square  root  of  the  input  signal  multiplied  by  the  square  root  of  10. 
Because  this  module  can  perform  three  different  functions,  it  is  very 
useful  for  simulating  certain  special  power  elements  such  as  induction 
motors  or  synchronous  machines. 

Comparator  Module 

This  module  consists  of  four  identical  comparators,  each  having  two 
inputs,  and  V^)  and  two  outputs  (C  and  C).  Because  the  signal  at 

the  input  terminal  serves  as  a reference  signal,  it  must  be  fixed 

at  a constant  potential.  The  other  input  signal,  fed  into  terminal  V^, 
can  be  varied.  The  sum  of  the  voltages  at  the  two  inputs  will  determine 
the  potential  states  of  the  outputs;  i.e.,  if  the  sum  of  these  two 
voltages  is  less  than  zero  the  output  C will  be  the  high  potential  state 
(or  5 V)  and  the  other  output  C will  be  the  low  potential  state  (or  0 
V).  If  the  sum  of  these  two  input  voltages  is  greater  than  zero,  the 
first  output  C will  be  the  low  potential  state  and  the  other  output  C 
will  be  the  high  potential  state.  This  change  in  the  potential  states 
of  the  outputs  can  be  used  to  control  an  analog  switch  to  connect  or 
disconnect  some  circuits  into  or  from  a simulated  power  system.  The  CTF 
of  the  simulated  power  system  will  vary  accordingly. 
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Function  Generator  Module 


This  module  is  designed  for  simulating  nonlinear  power  elements. 

The  method  used  in  this  simulation  is  similar  to  composing  a nonlinear 
curve  by  using  several  linear  segments.  If  the  x-coordinates  of  the 
starting  point  of  all  line  segments,  the  slopes  of  the  line  segments, 
and  the  y-coordinate  of  the  starting  point  of  the  first  line  segment  are 
known,  the  shape  of  this  curve  can  be  constructed.  The  function  generator 
module  is  designed  in  this  manner.  Figure  7 shows  that  the  module 
consists  of  10  sections  lined  in  sequence.  The  CTF  of  each  section  has 
been  programmed  in  such  a way  as  to  represent  a line  segment. 

The  main  components  of  this  module  are  a selecting  switch  SW;  a 
bias  potentiometer  ST;  a diode;  and  a slope  potentiometer  SL.  The 
selecting  switch  will  determine  the  sign  of  the  output  of  the  section  by 
switching  to  the  position  of  "+x"  or  "-x."  The  x-coordinate  of  a starting 
point  can  be  registered  into  a corresponding  section  by  setting  the 
input  voltage  equal  to  the  magnitude  of  x-coordinate;  i.e.,  by  adjusting 
ST  to  a position  at  which  the  diode  begins  to  conduct.  Once  the  diode 
is  conducting,  ST  and  SL  become  the  input  resistance  of  the  output 
amplifier.  Thus,  the  gain  of  the  amplifier  (or  the  slope  of  the  line 
segment)  is  controlled  by  SL.  Since  the  sections  in  the  module  are 
activated  in  sequence  as  the  value  of  x continuously  increases,  the  gain 
of  the  output  amplifier  is  a function  of  the  slopes  formed  by  the  first 
through  the  nth  line  segment;  thus,  the  registration  of  the  slopes  of 
those  line  segments  into  the  module  must  also  be  in  the  correct  sequence. 
The  registration  of  the  y-coordinate  of  the  starting  point  of  the  first 
line  segment  can  be  accomplished  by  setting  the  input  voltage  equal  to 
the  x-coordinate  of  the  starting  point  of  the  first  line  segment  and 
adjusting  the  voltage  of  an  additional  input,  y^^,  of  the  output  ampli- 
fier such  that  the  voltage  of  the  output  amplifier  is  equal  to  the  value 
of  the  y-coordinate. 
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SIMULATION  EXAMPLE 


; ! 

■ ! 


The  simulation  example  given  in  this  report  represents  the  actual 
power  distribution  system  at  the  Naval  Radio  Station  (NAVRADSTA(T)) , 

Isabela,  Puerto  Rico.  The  transients  at  this  station  have  been  studied 
and  reported  in  CEL  Technical  Note  TN-1239*  and  therefore  provide  a 
convenient  data  base  for  the  purpose  of  simulation.  The  power  distribu- 
tion system  shown  in  Figure  8 consists  of  one  substation  transformer  and 
three  distribution  transformers.  These  are  three-phase  transformers 
* with  A-connected  primaries  and  y-connected  secondaries.  The  name  plate 

information  of  these  transformers  are  shown  in  Figure  8.  Because  an  | 

underground  cable  is  used  for  transmitting  electrical  power  between  the 
substation  transformer  and  distribution  transformers,  the  circuit  diagram  | 

of  Figure  8 shows  capacitive  coupling  between  the  transmission  line  and  ! 

ground.  According  to  TN-1239,  the  value  of  these  capacitors  is  C = 0.58  pF.  ■ 

Total  load  of  this  power  system  is  400  kVA  at  an  0.8  power  factor. 

The  load  is  distributed  among  the  distribution  transformers  in  propor- 
tion to  their  ratings. 

The  electrical  transients  recorded  at  NAVRADSTA,  Isabela  were 
caused  by  lightning.  The  transient  wave  shape  is  shown  in  Figure  9,  and 
can  be  expressed  as 

-a, t -a,t  -o-t  -a„t 

aj^  t e + a2  e ® '’4  ^ (10) 

0.3079  ' 

0.7341 
0.2154 
5.925 
-19.062 
28.388 
-9.326 


*Naval  Civil  Engineering  Laboratory.  Technical  Note  N-1239:  Study  of 
electrical  power  distribution  system  transients  caused  by  lightning 
at  NAVRAOSTA(T) , Isabela,  Puerto  Rico,  by  J.  L.  Brooks  and  K.  Huang. 
Port  Hueneme,  Calif.,  Jun  1972. 


e,(t) 


where  a. 


= 


a.,  = 


a,  = 


This  transient  comes  from  the  three-phase  input  terminals  of  the 
substation  transformer,  thus,  the  transient  line-to-line  voltages  are 


i 

'i 


\ 


I 


"at 

-St  = 

(53)(2405)(6.8)  . ,,, 

35 

"bt 

-St  = 

(5)(2405)(6.8)  . 

55 

(11) 

®ct 

-St  = 

(48)(2405)(6.8)  _ 

55 

With  this  information,  the  power  system  can  be  simulated  as  follows. 
Sources  Simulation 

The  sources  which  excite  the  power  system  are  the  60-Hz  power 
source  and  the  transient  lightning  source.  These  sources  are  fed  into 
the  input  terminals  of  the  substation  transformer.  In  this  simulation 
these  sources  are  considered  as  infinitive  sources;  i.e.,  the  voltage 
waveform  of  the  sources  are  not  affected  by  the  impedance  of  the  power 
system.  With  this  consideration,  each  source  may  be  evaluated  individ- 
ually. 

(1)  Sinusoidal  Power  Source.  To  simulate  this  source,  a differen- 
tial equation  is  introduced. 

u)^  y = 0 (12) 

377 

of  this  simple  equation  is 

sin  U)t  + B cos  U)t  = C sin(u)t  + 6)  (13) 


where  U)  = 

The  solution 
y = A 


where  A aad  B are  arbitrary  constants  which  can  be  used  to  control  the 
initial  phase  angle  6 and  the  amplitude  C.  The  expressions  of  C and  6 
are  given  as 


(A^  + 

b2)1/2 

(lA) 

^ -1 
tan 

/B\ 

\aI 

(15) 

From  Equation  13  the  initial  conditions  of  y and  its  derivative  y are 

t = 0,  y = B and  y = A u)  (16) 

Based  on  Equations  12  through  16  the  sinusoidal  power  source  is  simulated 
as  shown  in  Section  1 of  Figure  10.  A scaling  factor  is  used  to 
prevent  overload  of  the  simulated  sinusoidal  source. 

(2)  Transient  Source.  The  Laplace  transformation  of  Equation  10  is 


E^(S) 


a.  a.  a.  a^ 

= + = + ■¥  — 

2 S + a S + a S + a 

(S  + Oj)  1 2 3 


Equation  17  consists  of  terms  in  the  form 


(17) 


E„(S) 


a 

S + a 


(18) 


Note  that  the  first  term  of  Equation  17  can  be  considered  to  be  tbe 
square  of  the  form  of  Equation  18.  The  time  domain  differential  equation 
of  Equation  18  is 


dt 


a e 


a 


= 0 


(19) 
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and  the  initial  condition  of  Equation  18  is 


t = 0,  = 0 (20) 

Therefore,  the  transient  waveform  expressed  in  the  equation  can  be 
simulated  by  programming  Equations  19  and  20.  The  analog  simulation  of 
these  equations  is  given  in  Figure  11.  In  a similar  fashion,  Equation  17 
can  be  programmed.  The  analog  simulation  is  given  in  Section  2 of 
Figure  10.  Scaling  factor  K2  is  used  to  prevent  the  analog  circuit 
Section  2 from  an  overload. 

(3)  Three-Phase  Source.  Since  the  high  voltage  side  of  the  substa- 
tion transformer  is  A-connected,  it  is  convenient  to  use  line-to-line 
voltage  (l*^,  and  (|»^  as  input  signals.  The  sum  of  these  three  signals, 
therefore,  must  equal  zero 

*•*3  '*’b  ***0  " ° 

where  6 = E,  - E 

^a  be 

0,  = E - E 

^b  c a 

<l»  = E - E, 

^c  a b 

The  voltages  E^,  E^,  and  E^  are  the  voltages  between  line  and  ground,  as 
shown  in  Figure  8.  The  60-Hz  source  is  a three-phase  source,  each  phase 
angularly  displaced  120  deg  from  the  other  two  phases.  Let  4*  be 

d 


By  use  of  Equations  21  through  23^  a three-phase,  60-Hz  voltage  can 
be  obtained  by  the  analog  schematic  diagram  shown  in  Section  3 of 
Figure  10. 

In  a similar  manner,  the  transient  signal  for  a three-phase  system 
can  be  simulated  by  employing  Equations  11  and  21.  The  analog  simulation 
of  this  is  shown  in  Section  3 of  Figure  10.  The  values  of  and  ^2 

o _ (48) (2405) (6. 8).  ^ - (5) (2405) (6. 8) 

Pi  ■ 33 * P2  ■ 33 

The  line-to-line  input  potential  of  the  substation  transformer  (ji^,  (j)^, 
and  represents  a potential  which  is  the  sum  of  the  AC  sinusoidal 
potential  and  the  transient  potential. 

Power  System 

The  main  task  in  the  simulation  of  a power  system  is  the  construc- 
tion of  a signal  flow  diagram  which  represents  the  power  system.  From 
this  diagram  an  actual  analog  circuit  can  be  connected.  Figure  12  is 
the  signal  flow  diagram  of  the  power  system  shown  in  Figure  8.  The 
input  signals  (J>^,  (|»^,  and  are  the  line-to-line  input  voltages  of  the 
substation  transformer.  These  input  signals  are  simulated  from  a source 
simulation  circuit  shown  in  Figure  10  and  are  stepped  down  by  the  ratio 
1/ng,  where 

ng  = 38,000/2,400  = 15.83 

Since  the  sources  are  assumed  to  be  infinitive  sources,  a feedback 

signal  to  the  simulated  source  circuit  is  not  required.  The  substation 

transformer,  therefore,  can  be  simply  represented  by  three  identical 

load  modules,  each  consisting  of  an  inductance  module  and  an  inductor 

series  resistor  module.  The  value  of  the  inductor  L can  be  calculated 

0 

from  the  nameplate  of  the  transformer 


15 


1 


‘o  ' 3 


(iTw)  (w)  = " 


The  resistor  R can  be  selected  such  that 
o 


R 

« 

U) 

U)  = 

0 

The  output  signals  of  these  modules  represent  the  output  currents  of  the 
transformer. 

The  next  power  element  on  the  signal  flow  diagram  will  be  three 
capacitors,  Y-connected  across  the  power  lines  to  the  ground.  Three 
capacitors  paralleled  with  resistor  modules  are  needed  to  simulate  the 
capacitors.  The  capacitance  of  these  modules  are 

C = 0.58  pF 
o 

and  resistance  of  the  modules  R are  chosen  such  that 


When  the  output  of  these  modules  represent  the  output  phase  voltages  of 
the  substation  transformer,  the  phase  voltages  must  be  converted  to  line 
voltages  through  three  summer  modules.  The  vector's  relationships 
between  phase  voltages  and  line  voltages  are  shown  in  Figure  13a.  These 
line  voltages  will  be  the  input  voltages  of  the  three-phase  transformers 
Tj^,  T2,  and  T^. 

Nine  transformer  modules  are  required  to  complete  the  simulation  of 
these  transformers.  Before  these  converted  line  voltage  signals  feed 
into  the  transformer  modules,  they  must  be  stepped  down  by  a ratio  1/n 
where 


n = 4,160/120  = 34.7 


I 


I 


1 

j 


I 

! 
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The  parameters  of  these  transformer  modules  are 


h = I5^  (3)  (^)  = 2 X 10-5  „ 

^2  ~ 500,000  ( 1,000 ) (377)  ■ 1.0  X 10  ^ H 

S ~ 225,000  (1,000)  (w)  " 1.6  X 10  ^ H 


and  L , , 
ml  ’ 

"m2’  "m3’  ®1’ 

and 

can  be  estimated  as 

"ml 

= 300 

"m2 

= 300 

"m3  = 300  L3 

«1 

_ 10 

10  P _ 10 

Lj  u)’ 

3 

CM 

3 1-3“* 

It  is  noted  in  Table  1 that  the  transformer  modules  have  two  input 
signals  (the  input  and  output  voltages  of  the  transformer)  and  two 
output  signals  (the  input  and  output  currents  of  ‘.he  transformer). 
Therefore,  these  transformer  modules  require  another  input  signal  to 
complete  input  information.  These  input  signals  for  the  nine  transformers 

^'^al'  ''a2’  '^a3’  '^bl’  ^2’  '^b3*  '^cl’  '^c2’  ^^03^  feedback  from  the 

outputs  of  the  respective  load  modules. 

One  of  the  output  signals  of  each  transformer  module  is  the  input 

current  of  the  simulated  transformer.  This  current  passes  through  the 

transformer  and  is  different  from  the  line  current.  Thus,  this  current 

must  be  stepped  down  by  a ratio  1/n  and  converted  to  line  current. 

Figure  13b  shows  the  relationship  between  the  currents  passing  through 

the  transformer  and  the  line  currents.  Line  currents  I',  IJ,  I’  will 

a’  b’  c 

feed  back  to  complete  the  current  information  of  the  whole  power  system. 
The  other  output  signal,  the  output  current  of  the  transformer,  passes 
through  load  modules  and  provides  the  output  voltage  for  each  phase  of 
these  transformers.  These  voltages  will  be  the  load's  response  to  the 
AC  sinusoidal  input  sources  as  well  as  to  any  transient  inputs. 
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Because  the  total  load  of  this  power  system  is  400  kVA  at  an  0.8 
power  factor  and  the  load  is  distributed  among  the  distribution  trans- 
formers in  proportion  to  their  ratings,  the  load  for  each  phase  can  be 
considered  as  a resistor  in  parallel  with  an  inductor.  Thus,  the  load 
connected  to  each  phase  of  these  transformers  can  be  calculated  as 
follows : 

(1)  The  load  on  each  phase  of  the  transformer  Tj^  is 


(400) (150) 

(150  + 500  + 225) (3) 


20,000 

875 


22.86  kVA 


The  current  in  each  phase  is 

, 22.86 


al 


120 

120 


= 190.5  A 

120 


al  0.8  I , ~ 0.8(190.5) 

al 


= 0.787  n 


uiL 


al  0.6  I 


120  _ 120 

0.6(190.5) 


= 1.05  fi 


r - 1-05 

al  ■ 377 


al 


= 2.785  mH 


and 


R I — R,  — R I 
al  bl  cl 


"al 


~ h)2  ~ h)3 


(2)  The  load  on  each  phase  of  the  transformer  T^  is 


(400) (500) 


(150  + 500  + 225) (3) 

The  current  in  each  phase  is 


I = = 635  A 

^a2  120  " 


= 76.19  kVA 
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0.8  I 


= 0.236  Q 


0.6  I 


= 0.315  ft 


= 0.835  mH 


®a2  ■ ^2  " ^c2 
^a2  " ^2  ” ^c2 

(3)  The  load  on  each  phase  of  the  transformer  is 


f 400) (225^ 


(150  + 500  + 225) (3) 

The  current  in  each  phase  is 


34,280  285 

120 


= 34.28  kVA 


0.8  I 


= 0.524  ft 


0.6  I 


= 0.699  ft 


'a3 

~ 377 

= 1 

a3 

■ ^3  ■ 

R - 
c3 

'a3 

” ^3  ” 

c3 

.855  mH 


(4)  All  chosen  values  of  resistances  in  the  load  module  of  an 
inductor  paralleled  with  a resistor  are  very  small,  compared  with  the 
reactance  of  the  corresponding  inductor. 
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After  all  parameters  in  the  signal  flow  diagram  are  defined,  the 
diagram  is  completed.  To  program  the  whole  power  system,  it  may  be 
necessary  to  scale  down  or  scale  up  the  signal  which  flows  in  the  simu- 
lated system.  The  scaling  factor  will  be  chosen  in  such  a way  that  the 
signal  magnitude  will  be  in  a range  between  1 and  9 V,  which  will  give 
maximum  accuracy. 


PROGRAMMING  TECHNIQUES 

The  basic  programming  techniques  for  PSS  are  the  same  as  the  program- 
ming techniques  for  the  analog  computer.  PSS  consists  of  many  modules 
which  have  been  preprogrammed  to  simulate  special  power  elements  such  as 
transformers  or  load  elements.  Thus,  to  avoid  unnecessary  progranming 
mistakes  and  confusion  on  the  interconnections  of  modules  and  the  adjust- 
ment of  their  performance  parameters , the  following  guidelines  provide 
aid  in  the  handling  of  these  modules. 

Circuit  Diagram 

A circuit  diagram  of  the  simulated  power  system  is  of  primary 
importance  for  accurate  simulation  because  all  connections  between  power 
element  modules  are  based  on  this  diagram.  The  diagram,  therefore, 
should  indicate  all  values  of  electrical  parameters  required  for  simula- 
tion of  the  power  system.  Some  circuit  diagrams  give  these  values  in 
ohms,  others  in  "pu"  (per  unit  impedance).  In  the  latter  cases,  these 
values  should  be  converted  to  resistance,  capacitance,  or  inductance 
because  this  conversion  will  make  all  parameters  independent  of  frequency. 
Having  completed  the  circuit  diagram,  the  next  step  is  to  construct  a 
signal  flow  diagram. 
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Signal  Flow  Diagram 


^ A power  element  has  at  least  one  pair  of  terminals  or  ports.  When  j 

the  element  is  energized,  two  parameters  can  be  measured  at  each  port: 

(1)  the  current  through  the  terminal  and  (2)  the  voltage  across  the  j 

terminals.  Theoretically,  either  one  of  these  parame*'“rs  can  be  desig-  j 

^ nated  as  the  input  signal  of  a power  element  module,  making  the  other  ! 

the  output  signal  of  the  module.  Once  a simulated  power  element  module  | 

■ designed,  its  input  and  output  signals  are  fixed,  and  the  other  power  j 

element  modules  connected  to  this  module  must  have  the  output  signal  of  ] 

, this  module  as  their  input  signal.  In  this  manner,  and  consistent  with 

the  circuit  diagram  of  the  simulated  power  system,  a signal  flow  diagram 
can  be  constructed.  In  the  example  in  this  report,  the  input  and  output 
signals  of  the  load  and  transformer  modules  were  preassigned  and  are 
listed  in  Table  1. 

I In  certain  cases,  however,  a difficulty  arises  in  constructing  a 

t signal  flow  diagram  in  the  prescribed  manner.  For  example,  consider  the 

i 

\ transformer  with  its  secondary  connected  to  a load  consisting  of  a 

I series  combination  of  an  inductor  and  resistor  as  shown  in  Figure  I4a.  j 

[ According  to  Table  1,  the  output  signal  of  the  secondary  of  the  trans-  | 

I former  module  is  output  current  i2,  while  the  input  signal  of  the  load  j 

' module  is  input  voltage  V.  Because  these  two  signals  are  different  in  | 

physical  units,  the  output  terminals  of  the  transformer  module  cannot  be  j 

I directly  connected  to  the  input  terminals  of  the  load  module.  To  over- 

! come  this  difficulty,  a dummy  element  must  be  inserted  between  these  two 

‘ elements.  The  dunmy  element  is  a resistor  of  relatively  high  resistance 

and  is  placed  across  the  secondary  of  the  transformer.  The  value  of 
resistance  must  be  so  chosen  as  to  have  only  a negligible  effect  on  the 
f simulated  circuit  when  connected.  Figure  14b  is  the  circuit  diagram 

■ , with  the  resistor  connected.  Figure  I4c,  then,  is  the  signal  flow 

1 diagram  of  the  circuit  in  Figure  I4b.  Thus,  the  dummy  element  serves 

! its  purpose-to  complete  the  signal  flow  diagram. 
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Scaling  Techniques 


Two  scaling  techniques  are  often  used  for  programming  a power 
system.  These  scaling  techniques  are  magnitude  scaling  and  frequency 
(time)  scaling.  All  simulating  modules  in  PSS  must  be  designed  with 
their  output  signals  in  the  range  of  1 to  10  V to  provide  the  most  accu- 
rate results.  Magnitude  scaling  techniques  will  scale  down  a large 
signal  to  the  proper  range  of  magnitude  by  use  of  a potentiometer  and  by 
scaling  up  a small  signal  by  adjustment  of  the  feedback  resistor  in 
the  operational  amplifier  of  the  module.  In  each  case  of  scaling  the 
signal  down  or  up,  the  scaling  factors  must  be  recorded  in  order  that 
the  real  magnitude  of  the  output  signal  can  be  determined. 

PSS  is  designed  to  operate  at  a center  frequency  of  60  Hz  with  a 
nominal  frequency  range  of  6 to  600  Hz.  The  simulation  of  high  frequency 
response  of  a power  system  beyond  600  Hz  will  result  in  decreased  accu- 
racy. In  order  to  overcome  this  defect  a frequency  scaling  technique 
can  be  employed.  The  technique  is  to  replace  the  high  frequency  source 
with  a source  whose  frequency  falls  within  the  nominal  frequency  range 
of  the  modules,  and  to  scale  all  values  of  inductances,  capacitances  and 
frequencies  as  indicated  by  Equations  24  and  25: 


U)  J 

L u)  = L ULj— ^ = (L  f„)  1 
0 0 o^  uij  o N 1 

1 _ 1 

2 

1 

1 II 

1 

s 

C u)  u)  ^ 1 

f f ' 

k C u) 

o 0 „ . 0 1 

^o  Wjj  0 N 1 

2 n f- 

1 

where  f^  = original  high  frequency 

fjj  = frequency  scaling  factor 

(24) 

(25) 


scaled  frequency  (between  6 and  600  Hz) 


“»o  = 

= 2nfjj 

u)  = 2nf 

^ = ^o/^N  = 

= original  inductance 

L = scaled  inductance 

= original  capacitance 

C = scaled  capacitance 


The  time  response  then  becomes 

^ ^0  (26) 

where  t^  = original  response  time 
t = scaled  response  time 

If  a linear  power  system  has  several  exciting  sources  whose  frequencies 
are  widely  varied,  the  superposition  theory  can  be  used  to  simulate  this 
system.  The  total  response  of  the  system  can  be  determined  by  summing 
the  responses  of  each  exciting  source  acting  alone.  Frequency  scaling 
techniques  may  be  required  to  let  the  exciting  frequency  fall  within  the 
nominal  frequency  range.  Then,  by  use  of  Equation  26  to  convert  all 
output  signal  waveforms  to  the  original  time  base,  the  total  response  in 
the  original  time  base  can  be  determined. 
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Characteristic  Transfer  Function  (CTF) 


To  simulate  a power  element  by  using  a simulating  module  of  PSS  as 
previously  mentioned,  it  is  necessary  to  match  the  CTF  of  the  module  to 
the  electrical  CTF  of  the  power  element.  As  indicated  by  Figure  A-1  and 
Equation  A-11  in  Appendix  A,  the  CTF  of  the  simulating  module  can  be 
expressed  as 


I*  = 


-A 


C S + 1 


S + 1 


(27) 


where  ({>  = CTF  of  the  simulating  module 

= function  of  the  wiper  position  on  potentiometer 

k^  = function  of  the  wiper  position  on  potentiometer  AR^ 


A = 


k. 

1 


B = kj  C 


(28) 

(29) 


If  the  electrical  CTF  of  the  power  element  to  be  simulated  is 


X = 


B'  S + 1 


(30) 


where  X is  the  electrical  CTF  of  the  power  element,  then  A'  and  B'  are 
known  constants. 

To  adjust  the  values  of  k^  and  k^  in  Equations  28  and  29,  the 
positions  of  the  wipers  of  potentiometers  AR^  and  AR^  are  adjusted  such 
that 


(31) 
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Then  the  simulating  module  is  analogically  identical  to  the  simulated 
power  element.  The  difference  in  sign  between  the  CTF  in  Equations  27 
and  30  is  not  significant  since  the  change  in  sign  of  (|i  can  easily  be 
accomplished  by  using  an  inverter. 

In  the  PSS,  however,  there  is  no  indication  of  the  values  of  k.  and 
k^;  therefore,  the  calibration  of  k^  and  k^  to  satisfy  Equations  31  and 
32  is  rather  difficult.  Fortunately,  this  difficulty  can  be  overcome  by 
checking  the  CTFs  in  Equations  27  and  30  with  VHF  signals  and  with  DC 
signals,  respectively.  For  VHF  signals. 


'•*|f»l  ~ ■ C k.  S 

A' 

* f»l  B'  S 
For  DC  signals. 


X 


f=0 


= A’ 


A 

B S 


(33) 

(34) 


(35) 

(36) 


From  Equations  33  through  36,  it  is  not  difficult  to  determine  a procedure 
for  matching  the  two  CTFs  by  using  the  procedure: 

(1)  Calculate  the  value  of  |x|  for  a fixed  high  frequency  f; 

i.e. , 


f»l 


A’  . 
B'  u)  ’ 


u»  s 2nf 


where  A' , B' , and  u»  are  known  values. 


(2)  Use  a unit  input  signal  at  the  fixed  high  frequency  to  excite 
the  simulating  module  and  adjust  such  that 

II  II  A' 

|*^|f»i  ~ |*|f»i  “ p~ur 

(3)  Calculate  the  value  of  ]x|j_q  at  DC  signal,  i.e., 


(4)  Use  a unit  input  DC  signal  to  excite  the  simulating  module  and 
adjust  such  that 

= |>'|f=0  ' *' 

(5)  Check  the  sign  of  i^,  if  (|)  and  x opposite  in  sign,  an  inverter 
is  needed  to  change  the  sign  of  i|i. 


CONCLUSIONS 

The  Power  System  Simulator  (PSS)  is  a powerful  tool  that  can  be 
used  to  represent  analogically  medium-sized  power  systems.  Once  pro- 
grammed, the  analog  computer  can  be  used  to  test  the  systems  for  tran- 
sients, load  response,  etc.,  by  \arying  the  parameters  of  the  modules. 

On  the  analog  simulation,  the  power  system  is  represented  analogi- 
cally by  various  hardware  system  modules.  Thus,  once  programmed  the 
latter  need  not  be  reprogrammed,  as  would  be  necessary  with  a software 
model.  Expansion  of  the  power  system  would  only  require  additional 
modules,  again  without  reprogramming  the  existing  analog  representation. 
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A definite  advantage  of  the  analog  simulation  is  that  voltages  and 
currents  can  be  recorded  with  meters  or  oscilloscopes  at  any  location  in 
the  system.  Unlike  a digital  computer,  where  each  parameter  must  be 
explicitly  specified  to  obtain  a measurement,  waveforms  can  be  recorded 
on  an  analog  simulation  by  moving  the  oscilloscope  probe  to  the  location 
in  question.  Thus,  signal  tracing  and  system  response  can  be  done  as 
easily  as  normal  troubleshooting  of  hardware  circuits. 

The  PSS  can  also  be  used  to  solve  other  engineering  problems. 

Basic  engineering  applications  of  this  equipment  are  described  in 
Appendix  B. 


27 


able  1.  List  of  Assigned  Input  and  Output  Signals  for  Programmable  Modules 


transformer  modules 


amplifier  modules 


voltage  reference  modules 


Figure  1.  Power  system  simulator. 


load  modules 


pot 


ntiometer  panels 


HIIICM  ClICII  III 


CX 


general  analog 
modules 


control  panel 
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Figure  2.  Simplified  circuit  diagram  of  amplifier  circuit. 


(a)  Resistor. 


(b)  Capacitor  parallel  (c)  Inductor  series 
with  a resistor.  resistor. 


Figure  3.  Typical  loads  in  power  system. 


Figure  5.  Signal  flow  diagram  of  a transformer  module. 


(c)  Square-root  device. 

Figure  6.  Configurations  of  the  interconnections  between  the  four 
^ terminals  of  the  multiplier. 
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(a)  Vector  Relationship  between  A-connection  transformer, 

line  voltage  and  phase  voltage. 


Figure  13.  Vector  and  current  relationships. 


(a)  Circuit  diagram  of  a transformer  connected  with  a load. 

Figure  14.  Circuit  diagram  of  a transformer  connected  with  a load  and  its 
signal  flow  diagram. 
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(c)  Signal  flow  diagram  of  (b). 


Appendix  A 


CTF  OF  SIMULATING  MODULES 

Figure  A-1  illustrates  the  basic  circuit  of  the  simulating  module. 

The  high  gain  operational  amplifier  in  the  circuit  provides  two  conditions: 


(1) 

and 


e.  « e.  and  e 
J 1 o 


(2)  i.  « i.  and  i- 
J 1 f 


For  simplification,  these  expressions  can  be  considered  as: 


(A-1) 

(A-2) 


Thus,  the  circuit  in  Figure  A-1  can  be  broken  into  two  separate  circuits 
as  shovm  in  Figure  A-2a  and  b.  The  loop  equations  of  the  circuit  in 
Figure  A-2a  can  be  written  in  matrix  form  as 


e. 

AR.  -X  AR. 

i 

1 

1 1 

a 

0 

-X  AR.  R.  + X AR. 

i . 

- 

11  1 

1 

X AR.  e. 

• 

X X 

X 

1 

AR.(R.  + X AR.)  - x^  AR.^ 

I X 1 

1 

i . 

^ 

1 1 

e. 

" R. 

K. 

1 

1 

X 

+ AR.(l-x)  ^ 

(A-3) 


(A-4) 
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where 


— + AR.(l-x) 

X 1 


(A-5) 


and  X is  proportional  to  the  wiper's  position  on  the  potentiometer 
its  value  is  between  zero  and  one. 

Since  R.  and  AR.  are  fixed,  value  K.  is  a function  of  x only. 
The  currents  in  the  circuit  of  Figure  A-2b  can  be  written  as 


^f  ~ ^fl  ^f2 


C S e 


^f2  ^f4  “ ^f3 


(A-6) 


(1-y)  AR^  i^3  + y ARj  i^^  = 


Rf  i^2  + y AR^  ij^  = 0 


Eliminating  i£^>  i£3  from  Equations  A-6  and  A-7  gives: 


(A-7) 


y aR£  i£^  % 

^f4  ■ ^f3  ■ ^f4  (1-y)  ARj  ~ (1-y)  AR^ 


1^,  e 

f4  o 

(1-y)  ■ (1-y)  ARj 

. ^f2  _ ^o 

" y(l-y)  AR,  (1-y)  AR 


«£  + y(i-y)  aRj 


+ (1-y)  ARj 


(A-8) 
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where 


Kf  = "Y  + (1-y)  AR^ 


(A-9) 


and  y is  proportional  to  the  wiper's  position  on  the  potentiometer  AR^; 
its  value  is  between  zero  and  one. 

Since  R^  and  AR^  are  fixed  value  is  a function  of  y only.  Then 


(A- 10) 


Combining  Equations  A-2,  A-4  and  A-10,  now  gives  the  CTF  of  simulating 
modules : 


-(Kj/K.) 

Kj  C S + 1 


-A 

B S + 1 


(A-11) 


APPENDIX  B , 

BASIC  ENGINEERING  APPLICATIONS  OF  PSS 


This  appendix  is  a reproduction  of  a manual  issued  by  the  Naval 
Civil  Engineering  Laboratory  on  the  application  of  the  equipment  dis- 
cussed in  this  report. 
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FOREWORD 


This  manual  is  arranged  for  the  purpose  of  providing  a framework  for  a 
training  course  on  the  engineering  applications  of  DC  750. 

Consequently,  the  emphasis  has  been  on  the  fundamentals  rather  than 
refinements. 

Furthermore,  since  the  readers  are  practicing  engineers  and  scientists, 
introductory  discussions  have  been  left  out,  or  reduced  to  bare  minimum. 
Complementary  discussions,  selected  based  on  the  general  interest  of  the 
participants,  will  be  presented  while  conducting  the  course. 
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I.  INTRODUCTION 


1.1  PHYSICAL  MODEL 

The  first  step  in  analyzing  an  engineering  problem  is  forming  a less  complex  and 
more  idealized  model  of  the  actual  physical  system;  the  outcome  of  this  step  is 
the  "physical  model." 

Example  1 — In  the  mechanical  system  shown  in  Figure  1 each  of  the  three 
components  has  mass,  elastic  deformation  and  internal  friction  (damping).  A 
"physical  model"  of  this  system  for  the  study  of  its  dynamic  behavior  is  formed 
by  assigning  to  each  element  its  predominant  quality.  Thus,  the  physical  model 
shown  in  Figure  2 is  obtained.  In  this  idealized  system  the  spring  and  the  dash- 
pot  have  no  mass,  there  is  no  elastic  deformation  in  the  dash-pot  and  the  mass, 
and  no  damping  anywhere  except  in  the  dash-pot. 


1.2  MATHEMATICAL  MODEL 

The  second  step  in  analyzing  an  engineering  problem  consists  of  expressing 
mathematically  the  cause-and-effect  relationships  of  each  component  and  the 
interactions  amongst  the  components.  The  "mathematical  model"  is  the  outcome 
of  this  step. 

Example  2 — The  physical  model  of  Example  1 results  in  the  following  mathematical 
model  (see  Figure  3 for  the  definition  of  symbols) : 
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Net  Force  on  M = MX 


F,  = KX 


F,  = BX 


Net  Force  on  M = F^  - (F,  + Fj) 


I ‘ 

IF2  l^i 


1*^2  “ *>=1 


IK  X B 


M = Mass 

K = Spring  Stiffness 

B = Dash-Pot  Friction 
Coefficient 

X = Displacement  of  M 


Figure  3 

Equations  (1  )-(3)  define  the  behavior  of  each  component,  and  Equation  (4) 
establishes  the  interaction  amongst  components. 

If  "Net  Force  on  M,"  F,  and  Fj  are  not  of  interest,  they  can  be  eliminated  and 
the  following  single  differential  equation  may  be  derived; 


MX  + BX-fKX  = F, 


1.3  COMPUTATION  STAGE 

The  third  step  in  analyzing  an  engineering  problem  is  evaluating  quantitatively  the 
response,  or  behavior  of  the  system,  under  a given  set  of  conditions;  this  is  the 
"computation  stage." 

Example  3 — In  the  mathematical  model  of  Example  2,  the  following  quantities 
are  given  in  a set  of  consistent  units: 


[ 

i 
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M = 1,000 


B = 6,000 

K = 144,000  , ^ 

I 

i 

At  an  instant  when  the  mass  M is  at  j 

X ^ X„  = 1 i 

and  has  an  upward  velocity  of  10,  i.e.,  ; 

1 

X ^ V„  = -10  i 

a downward  force  of  1 ,000,  i.e.,  I 

Fo  = 1,000,  ; i 

■ \ 

■ i 

is  applied.  How  does  the  mass  M move  from  this  instant  on?  j j 

The  purpose  of  "computation"  in  this  example  is  to  obtain  the  needed  quantity, 

X(t),from: 


Net  Force  on  M 
Fi 

F2 

Net  Force  on  M 


MX 
KX 
BX 
Fo  - 


(F,  + Fj) 


Mathematical  Model 


I 


(6) 


II.  METHODS  AND  MEANS  OF  COMPUTATION 


11.1  INTRODUCTION 

The  computation  stage  of  an  engineering  analysis  can  be  carried  out  in  a number  of 
ways.  The  choice  of  the  method  is  based  on: 

• Efficiency  of  computation, 

• Accuracy  requirement, 

• Mode  of  presentation;  numeric,  graphic,  etc., 

• Cost, 

• Availability  of  computer  and  conflict  with  other  users, 

• Frequency  of  need  to  compute  the  same  mathematical  model 
for  various  input  and  initial  conditions, 

• Number  of  parameters  and  their  range  of  variation, 

• Optimization  criteria  (when  applicable), 

• Whether  or  not  the  calculation  is  to  be  carried  out  within  a 
time  limit  while  events  are  actually  occurring  (e.g.,  events 
associated  with  controlling  a nuclear  reactor,  a satellite  or  a 
missile), 

• Technical  background  and  engineering  proficiency  of  the  analysis 
group,  etc. 

The  methods  and  means  of  fulfilling  the  "computation"  stage  of  an 
engineering  analysis  can  be  categorized  as  follows. 


II.2  DIRECT  SCALING  OF  PHYSICAL  SYSTEMS  (Prototype  Method) 

In  this  method  the  physical  system  is  scaled  dimensionally,  e.g.,  wind  tunnel 
evaluation  of  aircraft  structures. 

Example  4 — The  behavior  of  the  dynamic  system  described  by  the  mathematical 
model  of  Example  3 can  be  investigated  by  the  prototype  model  shown  in  Figure 
4 where: 


m = 

10- 

*M 

= 10-3 

X 1,000  = 1 

k = 

10- 

»K 

= 10-3 

X 144,000  = 144 

b = 

10- 

* B 

= 10-3 

X 6,000  = 6 

fo  = 

10- 

‘Fo 

= 10-3 

X 1,000  = 1 

Xo  = 

Xo 

= 

1 

Xo  = 

Xo 

= 

-10 

time  (prototype)  = time  (actual  system) 


It  can  be  easily  shown  that  the  mathematical  models  for  the  prototype  and  the 
original  system  are  identical,  the  initial  conditions  are  alike,  and  the  displacements 
are  equal  at  equal  elapsed  times. 


11.3  DIRECT  SIMULATION  OF  PHENOMENA  (Direct  Analog  Method) 

In  this  method  the  original  phenomena  is  replaced  by  phenomena  that  can  be  used 
more  conveniently,  and  implemented  with  less  cost  for  the  computation  purposes, 
e.g.,  photoelastic  stress  determination  of  complex  shapes. 

A.  Direct  Analogy  Using  Passive  Devices  (Passive  Analog  Method) 

In  this  approach,  simple,  elementary  devices  are  used  to  simulate  the  effects  of 
simple  elements  in  the  original  system,  e.g.,  an  electrical  capacitor  for  simulating 
a thermal  capacitance,  an  electric  resistor  for  a thermal  resistance,  etc. 

Example  5 — The  physical  model  considered  in  Example  2 can  be  investigated 
by  means  of  the  electrical  network  of  Figure  5.  In  fact,  the  mathematical  model 
for  this  network  is  (refer  to  Figure  5 for  definition  of  symbols) : 

* l/”'" 
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If  one  replaces  /vdt  by  a new  symbol,  say  x, 
the  above  expression  becomes; 

Cx  + (1/R)x  + (1/L)x  = Kt)  (10) 

Comparison  of  Equations  (5)  and  (10)  shows 
tfie  similarity  of  the  two  phenomena  and  the 
following  analogies  between  the  components 
of  the  two  systems; 


i 

Figure  5 


Mechanical  System 
Equation  (5) 


Electric  Model 
Equation  (10) 


M = Mass  *> 

B = Viscous  Friction 

K = Spring  Stiffness  o 

F = Force  -»• 

X = Velocity  o 


C = Capacity 

^ = Conductance 
R 

(l)  ~ Inverse  Inductance 

I = Current 
V = Potential 


B.  Direct  Analogy  Using  Active  Devices 

In  this  method  special  electronic  devices  with  amplifiers  and  other  active  circuits 
simulate  simple,  elementary  components  of  the  original  system.  The  main  limita- 
tion is  that  the  electronic  devices,  while  effective  and  trouble-free  for  integration, 
present  many  problems  in  differentiation;  furthermore,  these  devices  do  not 
operate  reversibly  as  their  counterparts  in  the  original  system  and  often  complex 
networks  are  needed  for  computer  modeling  of  the  reversibility. 

Example  6 — Arrange  a direct  analog  for  the  system  shown  in  Figure  6 with  the 
assumption  that  analog  devices  shown  in  Figure  7 are  available. 
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Let's  replace  the  active  devices  shown  for  (b),  (c)  and  (d)  by  the  simple  blocks 
shown  in  Figure  8;  label  the  blocks,  respectively,  R,  C-R  and  ■. 


V o — — J o-v 


Figure  8 


Using  these  blocks,  the  computer  model  shown  in  Figure  9 is  arranged. 


Figure  9 


As  it  can  be  seen,  this  is  an  "artificial"  direct  model,  one-to-one  correspondance  is 
lost  to  a great  extent,  and  modifications  in  the  physical  model  cannot  be  conven- 
iently implemented  in  the  computer  model. 


11.4  SOLVING  MATHEMATICAL  MODEL 

A.  Analytical  Solution  and  Numerical  Calculation  (Analytical  Method) 

In  this  method,  first,  the  mathematical  model  is  solved  analytically,  and  then 
evaluated  quantitatively.  Because  of  the  usually  complex  form  of  analytical 
.solutions  the  evaluation  is  almost  always  carried  out  numerically  (or  by  means 
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of  digital  computers).  In  a digital  computer,  functions  such  as  exponential, 
sine,  cosine  and  their  inverses,  square  root,  cube  root,  Bessel  functions,  etc., 
are  calculated  by  approximate  series  or  iterative  procedures  in  which  only 
principal  arithmetic  operations  are  employed. 

Example  7 — The  mathematical  model  of  Example  2 and  the  numerical  values 
of  Example  3 result  in  the  expression: 

X + 6X  + 144X  = 1 (Mathematical  Model)  (11) 


The  general  solution  of  the  above  equation  is: 


X(t)  = r:rr 


144 


+ X, 


1 e-3‘  sin(3  /iFt-  «/^) 

e-3»  sin{3  /ist-  >//) 

/i5 


where: 


+ X^  e-3*sin(3  yTTt) 

3/15 

4/  = tan{-  y 15) 


By  substituting  X and  X with  initial  values: 


Xo  = t, 


Xo  = -10, 


tor  any  given  t,  the  value  of  X(t)  is  numerically  calculated. 

Note  that  in  the  above  solution  time  is  a continuous  variable  and  X(t)  can  be 
evaluated  at  any  chosen  value  of  t. 

B.  Discrete  Digital  Computation  (Digital  Computer  Method) 

The  mathematical  models  of  engineering  systems  are  generally  integro-differential 
equations.  With  the  exception  of  very  simple  mathematical  models  where  ana- 
lytical solutions  are  readily  obtainable,  the  Discrete  Digital  Computation  is  used. 
In  this  method  the  derivatives  and  integrals  in  the  mathematical  model  are 
replaced  by  finite  difference  and  discrete  integration  formulae,  the  quantities 
of  interest  are  solved  for  and,  subsequently,  computed  digitally. 
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The  accuracy  of  calculated  results  depends  on  the  accuracy  of  finite  difference 
and  integration  formulae.  Instability  problems  resulting  from  calculations  are 
quite  frequent. 

Example  8 — The  mathematical  model  of  Example  7,  Equation  (1 1),  is 
considered  as  an  example.  By  using  the  finite  difference  formulae  (see 
Figure  10  for  the  definition  of  symbols): 

Xi  = ~ (Xm  - 2X;  + Xi,,) 

Xi  = ^ (Xi,i  - X;.,) 


Equation  (1 1)  becomes: 


Therefore,  using  Xj.,  and  Xj , one  can  compute  Xj^.,.  For  obtaining  X,, 
however,  a different  approach  is  to  be  used.  For  example,  since 
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— (X,  - X.,)  = - 10  (Given  initial  condition) 


X 


o 


one  obtains: 


X.1  = 20h  + X,  (14) 

and  Equation  (12),  written  for  i =0,  results  in: 

^ (X.1  - 2X„  + X,)  + 1 (X,  - X.i)  + 144X„  = 1 (15) 

which,  after  replacing  X^  with  the  initial  condition  1 and  X.,  with  Equation 
(14),  provides  X,,  i.e., 

- lOh  + 1 (16) 

Now  that  Xq  and  X-,  are  known.  Equation  (13)  can  be  used  to  compute 
Xj,  X3, ...,  successively. 

This  method  of  solution  suffers  particularly  from  discretization  errors  and 
limitations,  as  well  as  convergence  problems. 

C.  Continuous  Electronic  Computation  (Differential  Analyzer) 

In  this  method  use  is  made  of  electronic  devices  which  perform  the  principal 
arithmetic  operations  as  well  as  integration  with  respect  to  time.  Therefore, 
whenever  there  is  only  one  independent  variable  (usually  time)  in  the  system, 
the  outputs  of  the  system  can  be  computed  as  continuous  functions. 

This  method,  although  based  on  the  mathematical  model,  provides  a computer 
model  of  the  real  physical  system.  The  engineer  can  readily  associate  various 
signals  to  quantities  in  the  original  system  and  study  end-to-end  behavior  of 
the  system  while  the  computation  is  in  progress.  Furthermore,  rerunning  the 
computation  with  new  values  of  the  parameters,  or  local  modifications  of  the 
model  can  be  accomplished  readily  and  without  extensive  model  changes. 

Example  9 — Solve  the  mathematical  model  of  Example  7,  Equation  (11),  i.e., 

X + 6X  -t-  144X  = 1 (17) 
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with  the  initial  conditions: 


Xo  = 1,  Xo  = -10 

Since  X and  X are  successive  integrals  of  X,  the  graphical  representation  of 
Figure  1 1 always  holds. 


Figure  1 1 

However,  in  order  to  satisfy  the  mathematical  model.  Equation  (16),  one 

must  have  at  all  times;  * 

X = -6X  - 144X  + 1 (18) 

Thus,  X in  Figure  1 1 , instead  of  being  an  assumed  quantity,  must  be  generated  j 

from  X and  X as  shown  in  Figure  12.  , 
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In  a Differential  Analyzer,  electronic  devices  capable  of  integration  and 
summation  are  available;  however,  these  devices  introduce  a sign-change  also. 
The  Differential  Analyzer  model  shown  in  Figure  13  is  obtained  from 
Figure  12  directly  by  introducing  a sign  change  in  each  block  and  taking 
into  consideration  their  effects  throughout  the  network. 


III.  POWER  SYSTEM  SIMULATOR  DC  750 


III.1  CLASSIFICATION  OF  DC  750 

The  Power  System  Simulator  (DC  750)  is  a special-purpose  computer  designed 
primarily  for  the  simulation  of  power  systems. 

The  fundamental  component  of  DC  750  is  an  electronic  unit  called  "operational 
amplifier."  This  device  is  basically  a very  high  gain  amplifier  with  an  inherent 
sign  change. 

The  operational  amplifiers  in  DC  750  are  used  to  provide: 

a.  Summers  and  integrators. 

b.  Active  analog  modules  for  the  simulation  of  R,  CNR,  L-R  and  ideal 
transformers. 

With  the  present  components,  DC  750  is  capable  of  solving  engineering  problems 
of  moderate  size  by  the  two  methods; 

• Direct  analogy  by  means  of  active  devices  (see  Section  I I.2.B). 

• Electronic  differential  analyzer  computation  (see  Section  I I.3.C). 

The  Power  System  Simulator  provides  the  essential  equipment  and  circuits  (e.g., 
power  supplies,  control  switching  networks  and  instrumentation)  which  will  be 
required  for  active-passive  simulation  of  heat  transfer  and  dynamic  systems. 
Furthermore,  with  the  development  of  interface  and  logic  blocks,  complex 
problems  requiring  a hybrid  clfialog-digital  facility  can  be  solved. 


III.2  COMPONENTS  OF  DC  750 
A.  Power  Supplies 

Regulated  reference  DC  power  at  the  following  voltages: 

±10.  , ±1.  volts 

±100.  , ±10.  , ±1.  millivolts 

are  available.  Each  power  supply  is  protected  by  special  current  limiting  network 
and  overload  light  indicator.  All  power  supply  outputs  are  with  respect  to  a 
common  system  ground. 
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B.  Potentiometers 


There  are  two  potentiometer  panels,  each  panel  contains  16  ten-turn,  lOK-ohm 
potentiometers.  The  wiper  of  each  potentiometer  has  a 10  volt,  14  ma  current 
limiting  lamp  in  series  with  its  lead.  The  light  serves  as  an  overload  indicator 
as  well  as  a current  limiter  for  protecting  the  potentiometer. 


C.  Transfer  Trunks 

There  are  16  transfer  trunks  on  each  of  the  potentiometer  panels  numbered 
00-16.  The  even  numbered  trunks  have  outlets  on  the  control  panel;  all  of  the 
trunks  are  available  on  a distribution  panel  at  the  rear  of  the  computer  for 
external  connections  to  other  equipment. 

D.  Transistorized  Voltmeter 

A transistorized  voltmeter  is  available  on  the  control  panel.  A very  high  input 
impedance  amplifier  eliminates  the  loading  effect.  The  voltmeter  provides 
four  ranges;  30  V.,  10  V.,  3 V.  and  1 V. 


E.  Summers 

The  summers  are  formed  by  using  the  basic  operational  amplifiers  in  a feedback 
loop.  There  are  altogether  10  summers  available  in  5 dual  summer-integrator 
modules.  Figure  14  shows  input,  output  and  related  circuits  of  the  summers. 


F.  Integrators 

The  integrators  are  formed  by  using  the  basic,  chopper-stabilized,  operational 
amplifiers  in  a feedback  loop  with  provision  of  capacitor  in  the  feedback  branch. 
There  are  10  integrators,  one  in  each  of  the  dual  summer-integrator  modules. 
Figure  15  shows  input,  output  and  related  circuits  of  the  integrators. 


G.  General  Purpose  Module 

The  general  purpose  module  contains  three  units;  a summer,  an  inverter  and 
an  integrator  in  one  module.  Furthermore,  a number  of  feedback  resistor  and 
capacitor  selections  can  be  readily  made  by  means  of  the  range  switch  Rp/Ri  • 

Two  potentiometers  AR,  and  ARp  serve  to  adjust  the  integrator  input  and 
feedback  resistors  for  accurate  gains. 
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The  feedback  capacitors  are  C,  = 1 . /if  and  Cj  = 0.1  //f . The  range  switch 
utilizes  1 . /ifd  capacitor  for  the  right-half  of  the  dial  (switch  positions  marked 
10®,  10®,  10'*,  10^,  10^)  and  0.1  #i^d  capacitor  for  the  left-half  positions 
(marked  10®,  lO'*,  lO^,  10^,  10’). 

In  the  HG  (High  Gain)  position,  the  operational  amplifier  of  the  integrator  is 
essentially  in  open-loop  position. 

The  inputs,  outputs  and  related  circuits  are  shown  in  Figure  16. 

The  three  positions  of  the  panel  mode  switch  (toggle  switch  with  positions: 

R.  Only,  Norm.,  Cal)  and  the  range  switch  are  described  by  the  tables  in 
Figure  16  in  terms  of  the  switches  SW-1  and  SW-2  shown  in  the  figure. 

H.  Measurement  Network  and  Controls 

The  control  panel  contains  the  voltmeter  and  related  outputs,  reference  voltages, 
even-numbered  trunk  lines,  measurement  selection  matrix  switches,  mode 
control  switches  and  power  on-off  switch. 

The  function  of  measurement  selection  matrix  is  to  select  a desired  point  in 
the  computer  network  and  connect  it  to  the  READ  BUSS.  This  is  done  by 
pushing  the  push-button  switches  identifying  the  BAY  (1-3),  CAGE  (1-8), 
MODULE  (1-8)  and  the  particular  output  (A,  B,  C or  D)  in  the  module. 

The  mode-control  switches;  POT-SET,  RESET,  COMPUTE,  and  HOLD  are 
used  to  energize  the  switches  identified  PST,  RST  and  CMP  in  Figures  14  and 
1 5.  HOLD  switch  is  equivalent  to  having  neither  of  the  first  three  switches  on. 

The  power  ON-OFF  switch  is  a two-stop  push-button  switch;  the  second  stop 
is  spring-loaded  and  is  used  for  approximately  5 seconds  only  when  turning 
the  power  on.  To  turn  the  power  off,  the  switch  is  pushed  to  the  first  stop 
and  released.  The  green  light  of  "OFF"  sign  indicates  that  the  power  is 
available  to  the  computer.  The  green  light  must  be  on  before  attempting 
to  turn  on  the  computer. 


IV.  SCALING 


IV.1  INTRODUCTION 

In  solving  an  engineering  problem  on  DC  750,  all  variables  are  represented  by 
voltages  within  the  range  ±10  volts.  When  some  or  all  of  the  variables  in  the 
original  system  have  extremely  large  (or  extremely  small)  magnitudes,  scaling 
factors  are  required;  this  is  the  problem  of  "magnitude  scaling." 

Likewise,  when  the  dynamic  response  of  a system  is  to  be  investigated  on  DC  750, 
the  limited  frequency  response  of  the  computer,  as  well  as  the  desirable  time  dura- 
tion of  the  computer  runs,  necessitates  introducing  a time-scale;  this  is  the  problem 
of  "time-scaling." 


IV.2  MAGNITUDE  SCALING  (Normalizing  Method) 

Example  10 

[ Consider  the  problem  of  forming  a computer  model  for  the  mathematical  model: 

i 

X -I-  bx  + cx  = f(t)  (19) 

subject  to  the  initial  conditions 

i(0)  = V„,  x(0)  = x„ 

The  computer  model  is  shown  in  Figure  1 7 


Let  Ixi^x  > denote  estimates  of  maximum  magnitudes  of  x,  x 

and  X,  respectively,  and  be  the  maximum  of  the  given  function  f(t).  Then, 

the  new  variables 


X(t) 

- '<>  i?’ 

(21) 

X(t) 

■ '0 

(22) 

X(t) 

■ r!" 

(23) 

F(t) 

= 10 

lf(t)l,nax 

(24) 

are  functions  which  remain  within  ±10  volts  range.  , 

To  introduce  the  above  variables.  Equations  (19)  and  (20)  are  rewritten  as  follows: 


10  Xq 

X(0)  = (27) 

The  computer  model  for  Equations  (25)-(27)  is  shown  in  Figure  18 
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t i 


X(t)  = 10  ^ = 5 i(t) 


X(t)  = 10  ^ = x(t) 


and 


X(0)  = 20  x(0)  = 20  X (-0.5)  = -10. 
X(0)  = 5 x(0)  = 5x0  = 0 
Thus,  the  mathematical  model  becomes; 

X + 0.4X  + 0.8X  = 0 
The  computer  model  is  shown  in  Figure  19 


-X(0)  = 0 -X(0)  = -10 


(37 

(3£ 


(3£ 


Figure  19 


IV.3  TIME  SCALING 


Time  is  the  only  independent  variable  in  the  Differential  Analyzer.  Because  of  the 
frequency  response  limitations  of  the  computer,  or  undesirably  long  or  short  time 


duration  of  the  experiments  in  real  time,  scaling  of  time  is  often  necessitated. 
Time-scaling  is  introduced  at  either  mathematical  modeling  stage  or  computer 


A.  Time-Scaling  Mathematical  Model 


Example  12 

To  present  this  method  of  time-scaling,  consider  the  mathematical  model 

d^x  dx  , 

+ b—  -i-  cx  = f(t) 

dt2  dt 


with  the  initial  conditions; 


x{0)  = x„ 


x(0)  = V 

If  the  computer  time  is  denoted  by  t,  related  to  the  real  time  t by 

T = a t. 


then 


dx  _ ^ ^ ^ 

dt  dr  dt  dr 


dt2 


= a2 

dr  \ dt  / dt 


dr^ 


(40) 


(41) 

(42) 

(43) 

(44) 

(45) 


and  the  mathematical  model,  Equations  (40)-(42),  become: 


a' 


d^x 

dr^ 


+ ba 


dx 

dr 


-t-  c X = 


(46) 


x(0)  = Xo 


(47) 


dx 

dr 


r = 0 


1 ^ 

a dt 


t = o 


(48) 


71 


B.  Time-Scaling  Computer  Model 


The  only  "time-conscious"  component  in  DC  750  is  the  integrator.  The  time- 
constant  of  an  integrator  determines  its  rate  of  integration.  To  speed  up  the 
integration  process,  one  must  reduce  the  time  constant.  This  is  done  by 
reducing  the  value  of  either  the  input  resistor  or  the  feedback  capacitor. 

If  the  time-constant  of  the  integrator  is  reduced  by  reducing  the  input  resistor 
R without  changing  C,  since 


Integrator  Gain 


C 
R ' 


(49) 


the  gain  is  increased.  This  leads  to  an  alternative  way  of  implementing  time- 
scaling of  the  computer  model,  namely,  increasing  the  coefficient  factor  of 
input  potentiometers  by  the  required  factor,  or,  distributing  the  time  scale 
reduction  factor  between  the  gain  of  integrator  and  the  input  potentiometer 
coefficient. 


Example  13 

Consider  the  heating  problem  of  a well-insulated  body;  the  physical  model  is 
shown  in  Figure  20,  and  the  thermal  model  is  given  in  Figure  21 . 


Figure  20 

The  mathematical  model  is: 


-S'  = 


R 


1 


(50) 


where:  T = temperature  of  M 

C = thermal  capacity  of  M 


f 


R = thermal  resistance  of  insulation  + thermal 
resistance  of  boundary 

Tq  = initial  temperature  of  M 

Assume  the  following  numerical  values: 

C = 60Btu/°F 
R = 10°F.hr/Btu 
To  = -300° F 


Time-Scaling  of  the  computer  model  is  accomplished  by  changing  either  the 
integrator  time  constant,  or  the  integrator  input  potentiometer  coefficient,  or 
both. 
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models  shown  in  Figure  23  may  be  used. 


% 


r 


Assumptions: 


'1 


t 


1 — One-dimensional  flow  (no  corner  effects). 

2—  Neglect  ground  freezing  effect  on  properties. 

3—  Neglect  heating  of  cryogenic  due  to  cooling  of  the  ground 
under  the  tank. 

4—  Use  the  following  hypothetical  values: 

Tj  = Temperature  of  tank  (constant)  = -300°F 
= Temperature  of  underground  water  = 50°F 
T;,  = Atmospheric  temperature  = 70°F 
Kg  = 2 Btu/(hr.ft.°F) 

h^  = Convection  coefficient  to  underground  water  = 2.5 
Btu/(°F  ft2) 

hg  = Combined  convection  and  radiation  coefficient  to  the 
atmosphere  with  high  winds  = 3.2  Btu/(°F  ft^) 

Insulation  outside  the  tank  is  0.5  ft  thick  and  has  K|n  = 0.3 
Btu/(hr.ft.°F). 

Dimensions  are  shown  in  Figure  24. 

Thermal  Model — Figure  25 

G denotes  conductance,  i.e.,  inverse  of  resistance. 

Gi  = (2x1,256)/5  = 502.4 
G„  = 2.5x10x125.6  = 3,141 
Gg  = 3.2  X 1 ,256  = 4,020 
0.5  G,  = 251.2 

G,„  = (0.3  X 1,256)/0.5  = 753.6 
Mathematical  Model 


GaG, 

* it.-t,)  - 0 

0.5  G,  (T,  - Tj)  + 

Gw  <Tw  - h) 

(54) 

(T  -T) 
a -1-  G ' s 

+ 0.5Gi(T2 

- T,l  t G„(T„  - T,l  - 0 

(55) 

^ G +G, 

- T,) 

(56) 
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Let  Tj  = -300  be  represented  by -10  volts,  i.e., 


''ts  = X ^ = (-300)  X ^ = -10.  volts 
Then: 

' T.  X ^ - 70  X ^ - 2.333«,lts 
Vr.  - T„  X i - 60  X ^ - 1.667volts 

Vs  ” B IIT,  - T.)  G,G,/(G,  + G,ll  - 

Step  Three — Introduce  potentiometers  for  fractional  gains,  and  summer 
gains  for  factors  of  10.  See  Figure  28  for  the  result. 


Figure  28.  Final  Computer  Model  (Method  1) 
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Discussion: 


1—  Change  of  dimensions  and  parameters  without  change  of  model 
configuration. 

2—  Getting  additional  information  about  the  system  behavior. 

3—  A more  desirable  scaling  concept. 

4—  Reduction  of  network  size. 

Computer  Model — Method  2 (Passive  Analogy — 11.2-/ .) 

Step  One — Conversion  to  Electrical  Values. 

Let:  R = r(1/G) 

V = vT 


I = (v/r)(TG) 


The  units  are; 


[R]  = ohms,  [1/G]  = °F.hr/Btu 

[r]  = ohm/(°F.hr/Btu)  = 10.  Meg-ohm/(°F.hr/Btu) 

[V]  = volts,  [T]  = Of 

[v]  = volts/op  = — ^ (volt/OF) 

30Of  30 

[I]  = amp.,  [T  G]  = Btu/hr. 


[v/r]  = amp./{Btu/hr)  = pa/fBtu/hr) 

oUU 


The  electrical  counterparts  of  the  thermal  model  are: 


R,  = 10(1/502.4)  = 0.0199  Meg-Ohm  = 19.9K-Ohm 
= 10(1/3,141)  = 0.0031  B Meg-Ohm  = 3.18  K-Ohm 
R,  = 10(1/4,020)  = 0.00249  M = 2.49  K 
Ri„  = 10(1/753.6)  = 0.001328  M = 1.33K 


Step  Two — Electrical  Model,  Figure  29. 

2.333  volts 


Figure  29.  Computer  Model  (Method  2) 


B.  Transient  Heat  Flow 
Example  15 

Find  the  temperature-time  heating  curve  of  gasoline  in  a plastic  fuel  tank 
exposed  to  radiation  from  an  uncontrolled  fire,  Figure  30. 

The  thermal  system  of  Figure  30  may  be  replaced  by  the  physical  model  shown 
in  Figure  31 . 
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surrounding 


Simplified  Thermal  Model,  Figure  33 


Ro  = 


1 ^ ^ J_  / hr°F\ 

hpA  hpd)  hp  y Btu  j 


= 0.7,  Tp  = 150°F; 


Ro  = 


Ri  = 


Ri  = 


R2  = 


= 0.571  hr°F/Btu 


(1/16)/12  (1/2)/12 

0.1(1)  0.08(1) 


0.572  hr.OF/Btu 


Figure  33.  Thermal  Model 


+ TTiWi  = 10.52  hrOF/Btu 
0.08(1)  0.1(1) 


l/:j^V44)(0.5)  = 1.875  Btu/°F 


T,  = 


T,  = 


J_  ^ 1 

Rj  10.52 

Mathematical  Model 


0.095,  ^ 


= 0.534 


Qi  - (Tp  - T3)/Rp  = (Tp  - Tg)/R, 


(Tp  - TJ/R,  - (T.  - T,)/R2  = C — 


Tjt  = 0)  = T,  = 65° 


Measure  all  temperatures  with  respect  to  Ta  = Tf  = Tg(0)  = 65°,  then,  the 
mathematical  model  reduces  to: 

di  - Tp/Ro  = (Tp  - Tg)/R, 

(Tp  - Tg)/Ri  - Tg/R2  = C ^ 


Tg(t  = 0)  = 0 (Initial  Condition) 


Computer  Model — Method  1 (Diff.  An.,  II.3-C) 

Write  the  mathematical  model  as  follows: 


■ A " s;  ^)/if  " r) 


(63) 


■ U 


Ri 


iRi  ^2]  ^ 


dt 


(64) 


Tg(0)  = 0 


(65) 


J_  + JL  3.498 


= 0.286 


l/R, 

T 1 

R„  R 


1.748 

3.498 


= 0.5 


= 0.932.  + -^ 


C \R, 


C \ Ri  R, 


= 0.984 


Step  One — Assume  positive  gains  for  amplifiers.  (Figure  34) 

Tp  = 0.286  di  + 0.5  Tg 

t 

Tg  = f [0.932  Tp  - 0.984  Tg]  dt 


Tg(0)  = 0 
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Step  Three — Introduce  Scale-Factors,  Figure  36. 


LetTp  n^x  “ 500,  then; 


500  = v(10)  V = 50 


Degrees 


From; 


or 


Volts 

Tp  = 0.286  Qi  + 0.5  T, 

Tp  di  T 

— = 0.286  — -t-  0.5  — 

V V V 


one  concludes  that  qt  must  be  scaled  by  v also,  i.e.. 


J_ 

Comp.  Model  50 


Thermal 


6,700 

50 


= 134 


.f  i ^9) 


There  is  a problem:  the  highest  available  regulated  voltage  is  10.  volts. 
A partially  acceptable  solution  is  given  in  Figure  37. 
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Note:  134  x 0.286  = 38.3  = 10  x 0.383x10. 


Figure  37 


Computer  Model — Method  2 (Passive  Analogy,  II.2-A) 


Step  One — Thermal  Model,  Figure  38. 


Qi  = 
Ro  = 
Ri  = 

Rj  = 

C = 

"^g,  max 


6,700  Btu/hr 

0.571  hr°F/Btu 

0.572  hr°F/Btu 

10.52  hr^F/Btu 

1.875  Btu/°F 

1 = 220 

I expected 


"Ground"  represents  65°F 


Figure  38 


Step  Two — Select  scale  factors  and  form  Computer  Model,  Figure  39 


sw 


0,572M 


0.571  SM 


dummy 

load 

1 

> 

20  uf  < 

<1 1 

» < 

Then: 


Figure  39 

1 volt  20°F 
20/xfd  1.875  Btu/°F 
1.Meg<>  l.OF.hr/Btu 


1 volt 


= 1 fia  ~ 20  Btu/hr 


1 Meg 

20  nfd  X 1 Meg  = 20  sec  = 1 .875  hr 
or  1 sec  = 0.09375  hr  = 5.62  min 

(How  can  you  make  the  time  scale  1 sec  = 5 min?) 


Computer  Model  Components: 

' ■ 

Rq  = 0.571  x 1.  = 0.571  Megohms 
R,  = 0.572  Meg-ohms 
Rj  = 10.52  Meg-ohms 
C = 20/ifd 


Computation  is  accomplished  by  closing  the  switch  SW  and  recording 
Vg.  The  gas  temperature,  Tg,  is  obtained  by: 

Tg  = 20  Vg  + 55.  (OF) 


V.2  MECHANICAL  SYSTEMS 

A.  Impact  of  Falling  Bodies  (Drop  Test) 
Example  16 

The  package  shown  in  Figure  40  is 
dropped  to  impact  from  a height  h. 
Find  the  deceleration  of  M and  its 
displacement  with  respect  to  the 
container  after  the  impact.  The 
following  sets  of  parameters  are 
to  be  considered : 

M = 1 , 2gh  = 400 


a:  B = 4.5 

, 18 

, 22.5 

with  K = 81 

Figure  40 

b:  B = 6 

, 24 

CO 

o 

with  K = 144 

Assume  statical  deflection  negligible  compared  with  dynamical  deflection. 

Physical  Model,  Figure  41.  — ii — YX  ^ / 

At  the  impact  instant,  the  mass  M is  at  X = 0 ,,  J 

and  has  a velocity  = ||j_|  I 

//////  Y. 

Mathematical  Model  Figure  41 

MX  + BX  + KX  = 0,  X(0)  ^ V„  = = 20. 

or  X = -(B/M)X  -(K/M)X 

For  case  (a  1)  where;  B = 4.5,  K = 81,  we  have: 

X = -4.5X  -BIX,  B/M  = 4.5,  K/M  = 81. 
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Values  of  B/M  and  K/M  for  various  cases  are: 


Case; 

1 

2 

3 

■ 

B/M: 

4.5 

18 

22.5 

K/M: 

81 

81 

81 

B 

B/M: 

6 

24 

30 

K/M: 

144 

144 

144 

Computer  Model,  Figure  42. 


Vg  = 20  ic  • 0 


Note  that  (0)  = V^  = 20  volts  is  beyond  the  operating  voltage  range  of 
the  computer. 

The  maximum  deceleration  is  ~ 140  which  is  also  beyond  the  operating  voltage 
range  of  the  computer. 

To  solve  the  above  problems,  modify  the  computer  model  as  shown  in 
Figure  43. 
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Vic  - O’  Vo  - 2 V,c  - lOXo  “ 0 


B.  Shock  Absorber  Performance  Evaluation 
Example  17 

An  unbalanced  rotating  element  results  in  a disturbing  force  of  120  lb.  at  a 
repetition  rate  of  3,000  rpm.  The  machine  weighs  120  lb.  Four  vertical  springs 
with  stiffness  coefficients  1 ,500  Ib/in.  support  the  machine.  A single  shock 
absorber  with  a friction  coefficient  207  lb/(ft/sec)  is  proposed.  Find  the 
amplitude  of  vibration  and  maximum  resulting  acceleration. 

Physical  Model 
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Mathematical  Model 


MX  + 

X 

+ 

•X 

CO 

= f(t) 

(70) 

3.72  X + 207  X 

+ 72,000  X 

= 120sin314t 

(71) 

X + 56x  + 

19,300x  = 

32.3  sin  314 1 

(72) 

Frequency  of  unbalance  force  is  3,000/60  = 50  Hz  which  is  within  the 
nominal  frequency  range  of  DC  750  computer. 

Assume; 

1x1  max  = 5x10-^ 

= 0.2 

Ixl^,  = 100 


and  normalize  the  mathematical  model  as  follows: 

Ixl^^X  + 561x1^^  X + 19,300lxl,„xX  = 3.23 x 10  sin  314 1 (73) 

where; 


x/(5x10-*) 

= 2x103x 

(74) 

X = i/0.2 

= 5x 

(75) 

= x/100 

= 0.01  X 

(76) 

Rewrite  the  mathematical  model  as  follows. 


X = - 


56(0.2) 

100 


X - 


19,300(5)  10 
100 


1-4 


X + 


3.23 

100 


X 10  sin  314 1 


or 

X = -0.1 12  X - 0.0965  X + 0.0323  (10  sin  31 4t)  (77) 
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Computer  Model 


The  computer  model  is  shown  in  Figure  45. 

-sin  314t 


9 


* = ’00^  c-  « 

Figure  45 

Discussion: 

Computer  modeling  of  a simplified  version  of  VIBRA-LO. 


V.3  FLUID  SYSTEMS 

A.  Pipeline  Distribution  System 
Example  18 

In  the  pipeline  distribution  system  shown  in  Figure  46,  the  flow-rate,  q,  versus 
inlet-outlet  pressure  difference,  AP,  is  given  for  the  three  pipe  sections  in 
Figure  47.  The  pressure  at  the  pump  outlet,  P, , and  at  the  receiving  reservoirs, 
Pj  and  P3  are  known; 

P,  = SOpsi,  P2  = lOpsi,  P3  = 20psi 
Find  the  flow-rates  q2  and  q3. 
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Figure  47.  Pipe  Sections  Resistance  Characteristics 
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Draw  on  Figure  47  the  broken  lines  defined  by  the  algebraic  equations 

q = G AP  (78) 

G is  the  inverse  of  pipe  resistance. 

Physical  Model 

The  physical  model  is  shown  in  Figure  48. 


Mathematical  Model 

Let  P denote  the  pressure  at  the  branching  point  (Figure  48). 


The  flow-rate  equations  are. 

P,  - P = AP,,  G,  = f,(AP,).  q,  = GiAP,  (79) 

P - Pj  = APj,  Gj  = f2(AP2).  92  = GjAPj  (80) 

P - P3  = AP3.  G3  = f3(AP3).  q3  = G3AP3  (81) 

The  continuity  of  flow-rate  requires; 

q^  = q2  + q3  (82) 


Numerical  values  are:  P,  = 80,  Pj  = 10,  P3  - 20  and  the 
functional  relationships  given  in  Figure  47. 
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Computer  Model,  Method  1 (Diff.  An.,  II.3-C) 


Figure  49 


Since  the  magnitude  of  pressures  P,  and  P3  are  above  10,  scaling  of 
pressures  is  required.  Use  0.1  scale  factor,  i.e.,  (see  Figure  49). 


V, 

= 0.1  (P,) 

= 0.1(80) 

= 8 volts 

V2 

= 0.1  (Pj) 

= 0.1(10) 

= 1 volt 

V3 

= 0.1  {P3) 

= 0.1(20) 

= 2 volts 

This  scale  factor  will  affect  all  values  throughout  the  computer  model 
(Figure  50).  The  relationships  between  voltages  and  the  quantities 
they  represent  are 


Vi  = 

0.1 (P,)  = 

8 volts 

Vqi 

= 0.1  q. 

V2  = 

0.1  (P2)  = 

1 volt 

Vq2 

= 0.1  qj 

V3  = 

0.1  (P3)  = 

2 volts 

Vq3 

= 0.1  q3 

AV, 

= 0.1  AP,, 

AVj  = 0.1 

APj, 

t> 

< 

u 

It 

0.01  AP3 

Gi  = 

fi(APi)  = 

fidOAV,), 

G2  = 

fafAPj) 

= fjdOAVj), 

G3  = f3(100AV3) 
V = 0.1  P 
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Discussion:  Normalizing  method  of  scaling 


Effect  of  friction  reducing  agents 
Automation 


Computer  Model,  Method  2 (Passive  Analogy — 1 1.2-A) 

Step  1 — Scaling 

Lable  constant  G lines  in  Figure  47  by  corresponding  R = 1/G  values. 
From  these  values  of  R,  it  is  concluded  that  R^^g^  - ^5. 

Use  the  following  conversion  factors; 

r = 1,(X)0  ohm/1  (psi  sec/ft^)  (83) 

V = 0.1  volt/1  psi  (84) 

' “ ~ ~ afTip/(ft^/sec)  = 0.1  ma/(ft^/sec)  (85) 

r 1 ,uuu 

Step  2 — Forming  the  model 

Use  initial  guesses  for  pipeline  resistances  of  sections  1 , 2 and  3,  use  the 
above  conversions  and  form  the  computer  model  of  Figure  51 . 


Step  3 — Computer  Run 
a — Measure: 


AV,  = V,  - V 

AVj  = V - V2 

AV3  = V - V3 

b — Compute: 

AP,  = 10AV, 

APj  = IOAV2 
AP3  = IOAV3 

c — Use  Figure  51  and  AP,,  AP2,  AP3  to  determine  new  values  for 
R,,  R2  and  R3. 

d — Repeat  a,  b and  c until  R^ , R2  and  R3  match  values  read  from 
Figure  47  for  AP,,  AP2  and  AP3  (don't  forget  scale  factors). 

Discussion:  Automation 


1 

i 


B.  Dynamic  Response  of  a Pneumatic  Device 


Example  19 

The  physical  model  of  a component  in  a pneumatic  control  system  is  shown 
in  Figure  52.  Find  the  dynamic  response  of  this  device,  i.e.,  6(t)  for  a sudden 
change  of  the  line  pressure  from  to  P^  + P,.  The  numerical  values,  in 
lb,  slug,  in.,  sec  system,  are  as  follows: 

A = 4.  . L = 2.  , d = 4 , Po  = 1.375x10-® 

(l/R)  = 2x  10'®  (slug/sec)/(lb/in.^) 

K = 8 , P,  = 1.  , K is  the  stiffness  of  the  bellows. 
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Mathematical  Model 


For  small  x: 


* 


(P,  - P2)A  = Kx 


m,  = 


Pi  - P2 


Since 


mj  = A(L  - x)(Po  + P2) 


therefore  m2  s ALp2  - Apj,x 

From  ideal  gas  law: 

PVT'  = k,  7 = 1.4 

Since 

1 

P = V 

therefore 

P = kp’-^ 

From  the  above  expression: 


VI  .4 


Po  + P2 


klPo  + Pj)’"  = kp^  (1  ^ kp;"(l  + 1.4^ 


or 


'■  ■ {^) 


P2  ^3^2'  ^ I'^Pq/Po 


(86) 

(87) 

(88) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 
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For  air,  when  P is  in  Ib/in.^  and  p in  slug/in.^, 

K,  2 1.5x  10®  lb  in./slug  (95) 

at  sea  level  and  T = 60°F. 

Using  Equation  (94)  in  Equation  (89),  results  in: 

*2  = (AL/KJPj  - APoX  (96) 


and  Equation  (88)  becomes 

(AL/K3)P2  - Apoi  = (-^)(P,  - Pj) 

or,  when  rearranged : 


= L 


PoKa  . 1 K3 

" R 


Integration  of  both  sides  of  the  above  expression  results  in 
Po  Ka 


P2  = 


o Ka  Ka  r 

r"-  lit  /"■- 


(97) 


(98) 


(99) 


Computer  Model 

The  computer  model  is  arranged  by  using  Equations  (86),  (87)  and  (99), 
see  Figure  53. 

Numerical  values  are  as  follows: 


Ka 

RAL 


1.5x10®  X 2x10'® 
4x2 


0.375 


PoKa  1.375x10'®  X 1.5x10® 
L 2 


1.03 


1/d  = 0.25 

A scale  factor  of  1 volt/1  psi  may  be  used  for  converting  P, 
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NAVOCE  ANSYSCEN  Code  2010  San  Diego.  CA;  Code  3400  San  Diego  CA;  Code  409  (D.  G.  Moore).  San  Diego  CA; 
Ctxle  4473  Bayside  Library.  San  Diego.  C A;  Code  52  (H.  Talkington)  San  Diego  CA;  Code  5214  (H.  Wheeler).  San 
DiegoCA;  Code 5224 (R.Jones) San  DiegoCA;  Code 6505 (J.  Stachiw).  San  Diego.CA;  Code 6565 (Tech.  Lib.). 
San  Diego  CA;  Code  6700.  San  Diego.  C A;  Code  751 1 (PWO)San  Diego.  C A;  Code  81 1 San  Diego.  CA;  Research 
Lib..  San  Diego  CA;  SCE  (Code  6(500).  San  Diego  CA 
NAVORDSTA  PWO.  Louisville  KY 
NAVPETOFFCode  .30.  Alexandria  VA 
NAVPETRES  Director.  Washington  DC 

NAVPGSCOL  Code  1424  Monterey.  CA;  D.  Leipper.  Monterey  CA;  E.  Thornton.  Monterey  CA;  LCDR  K.C.  Kelley 
Monterey  CA 

NAVPHIBASECO.  ACB  2 Norfolk.  VA;  Code  S3T.  Norfolk  VA;  Harbor  Clearance  Unit  Two.  Little  Creek.  VA; 

OIC.  UCTONE  Norfolk.  Va 
NAVRADRECFAC  PWO.  Kami  Seya  Japan 

NAVREGMEDCEN  Chief  of  Police.  Camp  Pendleton  CA;  PWO  Newport  Rl;  PWO  Portsmouth.  VA;  SCE  (D.  Kaye); 

SCE  (LCDR  B.  E.  Thurston).  San  DiegoCA;  SCE.  Camp  Pendleton  CA;  SCE.  Guam 
NAVSCOLCECOFFC35  Port  Hueneme.CA;  CO.CodeC44A  Port  Hueneme. CA 

NAVSEASYSCOM  Code  0325.  Prtsgram  Mgr.  Washington.  DC;  Code  OOC  (LT  R.  MacDougal).  Washington  DC; 

Code  SEA  OOC  Washington.  DC 
NAVSEC  Code  6034  ( Library ).  Washington  DC 

NAVSECGRU  ACT  PWO.  Edzell  Scotland;  PWO.  Puerto  Rico;  PWO.  Torri  Sta.  Okinawa 
NAVSHIPREPFAC  SCE  Subic  Bay 

NAVSHIPYD;  CO  Marine  Barracks.  Norfolk.  Portsmouth  VA;  Code  202.4.  Long  Beach  CA;  Code  202.5  (Library) 
Puget  Sound.  Bremerton  WA;  Code  380.  (Woodroff)  Norfolk.  Portsmouth.  VA;  Code  400,  Puget  Sound;  Code 
4(X).()3  Long  Beach.  CA;  Code  404  (LT  J.  Riccio).  Norfolk.  Portsmouth  VA;  Code  410.  Mare  Is..  VallejoCA;  Code 
440  Portsmouth  NH;  Code  440.  Norfolk;  Code  440.  Puget  Sound.  Bremerton  WA;  Code  440.4.  Charleston  SC; 
Code  450.  Charleston  SC;  Code  453  (Util.  Supr).  VallejoCA;  L.D.  Vivian;  Library.  Portsmouth  NH;  PWD(Code 
4(K)).  Philadelphia  PA;  PW  D (LT  N.B.  Hall).  Long  Beach  CA;  PWO.  Mare  Is.;  PWO,  Puget  Sound;  SCE.  Pearl 
Harbor  HI; Tech  Library.  Vallejo,  C A 

NA  VSTA  CO  Naval  Station.  Mayporl  FL;  CO  Roosevelt  Roads  P.R.  Puerto  Rico;  Engr.  Dir..  Rota  Spain;  Maint. 
Cont.  Div..  Guantanamo  Bay  Cuba;  Maint.  Div.  Dir/Code  531.  Rodman  Canal  Zone;  PWD(LT  W.H.  Rigby). 
Guantanamo  Bay  Cuba;  PWO  Midway  Island;  PWO.  Keflavik  Iceland;  PWO,  Mayport  FL;  ROICC.  Rota  Spain; 
SCE.  Guam;  SCE.  San  Diego  CA;  SCE.  Subic  Bay.  R.P.;  Utilities  Engr  Off.  (LTJG  A.S.  Ritchie),  Rota  Spiiin 
NAVSTA  BISHOPS  K)INT  Harbor  Clear.  Unit  one.  Pearl  Harbor.  HI 

NAVSUBASE  ENS  S.  Dove,  Groton.  CT;  LTJG  D.W.  Peck,  Groton,  CT;  SCE.  Pearl  Harbor  HI 
NAVSUPPACT  CO.  Brooklyn  NY;  CO,  Seattle  WA;  Code  4.  12  Marine  Corps  Dist.  Treasure  Is.,  San  Francisco  CA; 

Code  413.  Seattle  WA;  LTJG  McGarrah.  VallejoCA;  Plan/Engr  Div..  Naples  Italy 
NAVSURFWPNCEN  PWO.  White  Oak.  Silver  Spring,  MD 
NAVTECHTRACEN  SCE.  Pensacola  FL 
NAVTORPSTA  Keyport.  WA 

N AVWPNCEN  C(xle  26.36  ( W.  Bonner).  China  Lake  CA;  PWO  (Code  26).  China  Lake  CA;  ROICC  (Code  702).  China 
Lake  C A 

NAVWPNEVALFACTechnical  Library.  Albuquerque  NM 

NAVWPNSTA  (CTebak)Colts  Neck.  NJ;  Code  092.  Colts  Neck  NJ;  Code  092A  (C.  Fredericks)  Seal  Beach  CA;  ENS 
G.A.  Lowry,  Fallbr<«>k  CA;  Maint.  Control  Dir.,  Yorktown  VA;  PW  Office  (Code09CI ) Yorktown.  VA;  PWO. 
Seal  Beach  CA 

NAVWPNSUPPCEN  Code  09  (Boennighausen)  Crane  IN 
NAVXDIVINGU  LT  A M.  Paris!.  Panama  City  FL 
NCBU  405  OIC.  San  Diego.  CA 
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NORDA  Code  410  Bay  Si  Louis,  MS 

NCBCCEl.  (CAPT  N.  W.  Petersen),  Port  Hueneme,  CA;  CEL  AOIC  Port  Hueneme  CA;  Ccxle  10  Davisville,  Rl; 

Code  15.',  Port  Hueneme  CA;  Code  1.56,  Port  Hueneme,CA;  Ctsde  251 1 1 Port  Hueneme,  C A;  Code  400,  Gulfport 
MS;  NESOCode251  P.R,  Winter  Port  Hueneme,  CA;  PW  Engrg,  Gulfport  MS;  PWOfCode 80) Port  Hueneme, 
CA,  PWO,  Davisville  Rl 
NCBU411  QIC,  Norfolk  V A 
NCR  20,  Commander 
NCSO  BAHRAIN  Security  Offr,  Bahrain 

NMCB  133(ENST.W.  Nielsen);  5,  Operations  Dept  .;  Forty,  CO;  THREE,  Operations  Off. 

NOAA  Libraries  Div.  - D823.  Silver  Spring,  MD 
NORDA  Code  440  (Ocean  Rsch,  Off)  Bay  St.  Louis.  MS 

NRL  Code  8400(3.  Walsh),  Washington  DC;  Code  8441  (R.A.  Skop).  Washington  DC;  Rosenthal.  Code  8440.  Wash. 
DC 

NSC  Code  54.1  (Wynne),  Norfolk  VA 
NSDSCE.  Subic  Bay.  R.P. 

NTC  Code  54  (ENS  P.  G.  Jackel),  Orlando  FL;  Commander  Orlando.  FL;  SCE  Great  Lakes.  IL 
NAVOCEANSYSCEN  Hawaii  Lab(D.  Moore).  Hawaii 

NUSCCode  131  New  London.  CT;  Code  EAI23(R.S.  Munn).  New  London  CT;  Ctxle  S332.  B-80  (3.  Wilcox);  Code 
SB  331  (Brown).  Newport  Rl;  Code  TA131  (G.  De  la  Cruz).  New  London  CT 
OCEANAV  Mangmt  Info  Div.,  Arlington  VA 
OCEANSYSLANT  LT  A.R.  Giancola,  Norfolk  VA 

OFFICE  SECRETARY  OF  DEFENSE  OASD(l&L)  Pentagon  (T.  Casberg).  Washington  DC 

ONRCDR  Harlett.  Boston  MA 

NORDA  Code  440(Ocean  Rsch.  off)  Bay  St.  Louis.  Ms 

ONR;  BROFF.  CO  Boston  M A;  Code  221 . Arlington  VA;  Code  481 . Arlington  VA;  Code  700F  Arlington  VA;  Dr.  A. 

Laufer.  Pastidena  CA 
PACMISRANFACCO.  Kekaha  HI 

PMTC  Code  4253-3,  Point  Mugu,  CA;  Pat.  Counsel.  Point  Mugu  CA 

PWC  ENS  3.E.  Surash.  Pearl  Harbor  HI;  ACE  Office  (LT3G  St.  Germain)  Norfolk  VA;  CO  Norfolk.  VA;  CO.  Great 
Lakes  IL;  Code  1 16(LT3G.  A.  Eckhart)Great  Lakes,  IL;  Code  120.  Oakland  CA;  Code  1 20C  ( Library ) San  Diego. 
CA;  Code  128.  Guam;  Code  200,  Great  Lakes  IL;  Code  200.  Oakland  CA;  Code  220  Oakland.  CA;  Code  220.1. 
Norfolk  VA;  Code  .30C (Boettcher) San  Diego. CA;  Code40(C.  Kolton) Pensacola.  FL;  Code42B(R.  Pascuii). 
Pearl  Harbor  HI;  Ctxle  .'05  A (H . Wheeler);  Code  680.  San  Diego  CA;  OIC  CBU-405.  San  Diego  C A;  XO  Oakland. 
CA 

spec  Code  1 22B.  Mechanicsburg.  PA;  PWO  (Code  1 20)  Mechanicsburg  PA 
TVA  Smelser.  Knoxville.  Tenn. 

U.S.  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

US  GEOLOGICAL  SURVEY  Off.  Marine  Geology.  Piteleki.  Reston  VA 

USAF  Maj.  Riffel,  Rumstein.  Germany 

USAF  REGIONAL  HOSPITAL  Fairchild  AFB.  WA 

USCG(G-ECV)  Washington  Dc;  (G-ECV/61)(Burkhart)  Washington.  DC;  (G-MP-3/USP/82)  Washington  Dc; 

G-EOE-4/61  (T.  Dowd),  Washington  DC 
USCG  ACADEMY  LT  N.  Siramandi.  New  London  CT 

USCG  R&D  CENTER  COGroton.  CT;  D.  Motherway.  Groton  CT;  LT3G  R.  Dair,  Groton  CT;  Tech.  Dir.  Groton.  CT 
USNA  Ch.  Mech.  Engr.  Dept  Annapolis  MD;  Energy-Environ  StudyGrp.  Annapolis.  MD;  Engr.  Div.  (C.  Wu) 
Annaplolis  MD;  Environ.  Prot.  R&D  Prog.  (3.  Williams).  Annapolis  MD;  Ocean  Sys,  Eng  Dept  (Dr.  Monney) 
Annapolis.  MD;  PWD  Engr.  Div.  (C.  Bradford)  Annapolis  MD 
AMERICAN  UNIVERSITY  Washington  DC(M.  Norton) 

ARIZONA  Suite  Energy  Programs  Off..  Phoenix  AZ 

BONNEVILLE  POWER  ADMIN  Portland  OR  (Energy  Consrv.  Off..  D.  Davey) 

CALIF.  MARITIME  ACADEMY  Vallejo.  CA  (Library) 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY  Pasadena  CA  (Keck  Ref.  Rm) 

CALIFORNIA  STATE  UNIVERSITY  LONG  BEACH.  CA  (CHELAPATI) 

CATHOLIC  UNIV.  Mech  Engr  Dept.  Prof.  Niedzwecki.  Wash..  DC 
CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  & M(X)RE  LIBRARY  LOS  ANGELES.  CA 

DUKE  UNIV  MEDICAL  CENTER  B.  Muga.  Durham  NC 

FLORIDA  ATLANTIC  UNIVERSITY  BOCA  RATON.  FL(MC  ALLISTER) 
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FLORIDA  ATLANTIC  UNIVERSITY  Bwa  Raton  FL  (W.  Tessin) 

FLORIDA  TECHNOLOGICAL  UNIVERSITY  ORLANDO.  FL  (HARTMAN) 

FOREST  INST.  FOR  (X'EAN  & MOUNTAIN  Carson  City  NV  (Studies  - Library) 

FUEL  & ENERGY  OFFICE  CHARLESTON.  WV 

HAWAII  STATE  DEPT  OF  PLAN.  & ECON  DEV.  Honolulu  HI  (Tech  Info  Ctr) 

ILLINOIS  STATE  GEO.  SURVEY  Urbana  11. 

INDIANA  ENERGY  OFFICE  Energy  Group.  Indianapolis.  IN 
VIRGIN  lA  INST.  OF  MARINE  SCI.  Gloucester  Point  VA  (Library) 

KEENE  STATE  COLLEGE  Keene  NH  (Cunningham) 

LEHIGH  UNIVERSITY  BETHLEHEM.  PA  (MARINE  GEOTECHNICAL  LAB..  RICHARDS);  Bethlehem  PA 
(Linderman  Lib.  No.30.  Flecksteiner) 

LIBRARY  OF  CONGRESS  WASHINGTON.  DC  (SCIENCES  & TECH  DIV) 

LOUISIANA  DIV  NATURAL  RESOURCES  & ENERGY  Dept,  of  Conservation.  Baton  Rouge  LA 
LOW  COUNTRY  REG.  PLAN.  COUNCIL  YEMASSEE.  SC(BAGGS) 

MAINE  MARITIME  ACADEMY  (Wyman) Castine  ME;  CASTINE.  ME  (LIBRARY) 

MAINE  OFFICE  OF  ENERGY  RESOURCES  Augusta,  ME 
MISSOURI  ENERGY  AGENCY  Jefferson  City  MO 

MIT Cambridge  MA;  Cambridge  M A (Rm  10-500.  Tech.  Reports.  Engr.  Lib.);  Cambridge.  MA  (Harleman) 

MONTANA  ENERGY  OFFICE  Helena.  MT 

NATL  ACADEMY  OF  ENG.  ALEXANDRIA.  VAISEARLE,  JR.) 

NEW  MEXICO  SOLAR  ENERGY  INST.  Dr.  Zwibel  Las  Cruces  NM 
NY  CITY  COMMUNITY  COLLEGE  BROOKLYN.  NY  (LIBRARY) 

NYS  ENERGY  OFFICE  Library,  Albany  NY 
OREGON  Salem  OR  (Suite  Office  of  Allocation  & Conserv.) 

OREGON  STATE  UNIVERSITY  (CE  Dept  Grace)  Corvallis.  OR;  Corvalis  OR  (School  of  Oceanography) 
PENNSYLVANIA  STATE  UNIVERSITY  STATE  COLLEGE.  PA  (SNYDER);  State  College  PA  (Applied  Rsch  Lab) 
POLLUTION  ABATEMENT  ASSOC.  Graham.  Corte  Madera.  CA 
PURDUE  UNIVERSITY  Lafayette.  IN  (CE  Engr.  Lib) 

CONNECTICUT  Hanford  CT(Dept  of  Plan.  & Energy  Policy) 

SCRIPPS  INSTITUTE  OF  OCEANOGRAPHY  LA  JOLLA.  CA  (ADAMS);  San  Diego.  CA  (Marina  Phy.  Lab.  Spiess) 
SEATTLE  U Prof  Schwaegler  Seattle  WA 

SOUTHWESTRSCH  JNSTJ.  Maison.  San  Antonio  TX,  R.  DeHart.  San  Antonio  TX 
STANFORD  UNIVERSITY  Engr  Lib.  Stanford  CA;  STANFORD. CA  (DOUGLAS) 

STATE  HOUSE  AUGUSTA.  ME  (MAINE  STATE  FUEL  ALLOC  & CONSERV.  OFF.) 

STATE  UNIV.  OF  NEW  YORK  Fort  Schuyler.  NY  (Longobardi) 

TEXAS  A&M  UNIVERSITY  COLLEGE  STATION.  TX  (CE  DEPT);  College  Station  TX  (CE  Dept.  Herbich) 
UNIVERSITY  OF  CALIFORNIA  BERKELEY.  CA  (CE  DEPT.  GERWICK);  BERKELEY.  CA  (OFF.  BUS.  AND 
FINANCE,  SAUNDERS);  Berkeley  CA(E.  Pearson);  LIVERMORE.  CA  (LAWRENCE  LIVERMORE  LAB. 
TOKARZ);  La  Jolla  CA  (Acq.  Dept.  Lib.  C-075A);  Off-Campus  Facs  & Energy  Conserv.,  Davis,  CA;  Off-Campus 
Facs  & Energy  Conserv..  Davis.  CA;  SAN  DIEGO.  CA.  LA  JOLLA,  CA  (SEROCKI);  Vice  President.  Berkeley, 
CA 

UNIVERSITY  OF  DELAWARE  Newark,  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  Dr  Chiu  Honolulu.  HI;  HONOLULU.  HI  (SCIENCE  ANDTECH.  DIV.) 

UNIVERSITY  OF  ILLINOIS  Metz  Ref  Rm.  Urbana  IL;  URBANA.  IL  (DAVISSON);  URBANA.  IL  (LIBRARY); 
URBANA.  IL(NEWARK) 

UNIVERSITY  OF  KANSAS  Kansas  Geological  Survey.  Lawrence  KS 
UNIVERSITY  OF  MASSACHUSETTS  (Heronemus).  Amherst  MA  CE  Dept 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  MI  (Richart) 

UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE  (Ross  Ice  Shelf  Proj.) 

UNIVERSITY  OF  NEW  HAMPSHIRE  DURHAM.  NH  (LAVOIE) 

UNIVERSITY  OF  RHODE  ISLAND  KINGSTON.  RI  (PAZIS);  Narragansett  Rl  (Pell  Marine  Sci.  Lib.) 

UNIVERSITY  OF  SO.  CALIFORNIA  Univ  So.  Calif 

UNIVERSITY  OF  TEXAS  Inst.  Marine  Sci  ( Library ).  Port  Arkansas  TX 

UNIVERSITY  OF  TEXAS  AT  AUSTIN  AUSTIN.  TX  (THOMPSON) 

UNIVERSITY  OF  WASHINGTON  Dept  of  Civil  Engr  (Dr.  Matlock).  Seattle  WA;  SEATTLE.  WA  (APPLIED 
PHYSICS  LAB);  SEATTLE.  WA  (OCEAN  ENG  RSCH  LAB.  GRAY);  SEATTLE.  WA  (PACIFIC  MARINE 
ENVIRON.  LAB..  HALPERN);  Seattle  WA  (E.  Linger);  Seattle.  W A Transportation.  Construction  & Geom.  Div 
UNIVERSITY  OF  WISCONSIN  Milwaukee  Wl  (Ctr  of  Great  Lakes  Studies) 
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URS  RESEARCH  CO.  LIBRARY  SAN  MATEO.  CA 

VENTURA  COUNTY  ENVIRON  RESOURCE  AGENCY  Ventura.  CA  (Melvin) 

VERMONT  STATE  ENERGY  OFFICE  MONTEPELIER.  VT(DIRECTOR) 

AGBABIAN  ASSCK.  C.  Bugge.  El  Segundo  CA 

AMETEK  Offshwe  Res.  & Engr  Div 

AMSCO  Dr.  R.  McCoy.  Erie.  PA 

ARCAIR  CO.  D.  Young.  Lancaster  OH 

ARVID  GRANT  OLYMPIA.  WA 

ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMITH) 

USNA  Sys.  Engr.  Dept  (R.  McCoy)  Annapolis.  MD 
BECHTEL  CORP.  SAN  FRANCISCO.  CA  (PHELPS) 

BELGIUM  HAECON.  N.V  .Gent 

BRITISH  EMBASSY  Sci.  & Tech.  Dept.  (J.  McAuley).  Washington  DC 
BROW  N & R(X)T  Houston  TX  (D.  Ward) 

CANADA  Adrian,  ^nderson  & Assoc..  Winnipeg;  Can-Dive  Services(English)  North  Vancouver;  Library.  Calgary. 
Alberta;  Lockheed  Petro.  Serv.  Ltd.  New  Westminster  B.C.;  Lockheed  Petrol.  Srv.  Ltd..  New  Westminster  BC; 
Nova  Scotia  Rsch  Found.  Corp.  Dartmouth.  Nova  Scotia;  Surveyor.  Nenninger  & Chenevert  Inc..  Montreal 
CAYWOOD.  NOPP.  W ARD.  AlA  Sacramento.  CA 
CHEMED  CORP  Lake  Zurich  IL  (Dearborn  Chem.  Div. Lib.) 

CHEVRON  OIL  FIELD  RESEARCH  CO.  LA  HABRA.  CA  (BROOKS) 

COLUMBIA  GULF  TRANSMISSION  CO.  HOUSTON.  TX  (ENG.  LIB.) 

CONCRETE  TECHNOLOGY  CORP.  TACOMA.  WA  (ANDERSON ) 

DESIGN  SERVICES  Beck.  Ventura.  CA 
DILLINGHAM  PRECAST  F.  McHale.  Honolulu  HI 
DIXIE  DIVING  CENTER  Decatur.  GA 

DRAVO  CORP  Pittsburgh  PA  (Giannino);  Pittsburgh  PA  (Wright) 

DURLACH.  O'NEAL.  JENKINS  & ASSOC.  Columbia  SC  • 

NORWAY  DET  NORSKE  VERITAS  (Library ).  Oslo 
EVALUATION  ASSOC.  INC  KING  OF  PRUSSIA.  PA(FEDELE) 

EXXON  PRODUCTION  RESEARCH  CO  Houston  TX(  A.  Butler  Jr) 

FORD.  BACON  & DAVIS.  INC.  New  York  (Library) 

FRANCE  Roger  LaCroix.  Paris 

GLIDDEN  CO.  STRONGSVILLE.  OH  (RSCH  LIB) 

GLOBAL  M ARINE  DEVELOPMENT  NEWPORT  BEACH.  CA  (HOLLETT) 

GOULD  INC.  Shiidy  Side  MD(Ches.  Inst.  Div..  W.  Paul) 

GRUMMAN  AEROSPACE  CORP.  Bethpage  NY  (Tech.  Info.  Ctr) 

HUGHES  AIRCRAFT  Culver  City  CA(Tech.  Doc.  Ctr) 

MAKAl  OCEAN  ENGRNG  INC.  Kailua.  HI 
JAMES  CO.  R.  Girdley.  Orlando  FL 

LAMONT-DOHERTY  GEOLOGICAL  OBSERV.  Palisades  NY  (Selwyn) 

LOCKHEED  MISSILES  & SPACE  CO.  INC.  L.  Trimble.  Sunnyvale  CA;  Mgr  Naval  ArchiS:  Mar  Eng  Sunnyvale. 

CA;  Sunnyvale  CA  (Rynewicz);  Sunnyvale.  CA  (Phillips) 

MARATHON  OIL  CO  Houston  TX  (C.  Seay) 

MCDONNEL  AIRCRAFT  CO.  Dept  .“iOl  (R.H.  Fayman).  St  Louis  MO 

MEDERMOTT&CO.  Diving  Division.  Harvey.  LA 

MOBIL  PIPE  LINE  CO.  DALLAS.  TX  MGR  OF  ENGR  (NOACK) 

NEWPORT  NEWSSHIPBLDG  & DRYDOCKCO.  Newptirt  News  VA  (Tech.  Lib.) 

NORWAY  A.  Torum.  Trondheim;  DET  NORSKE  VERITAS  (Roren)  Oslo;  1.  Foss.  Oslo 
(XTEAN  DATA  SYSTEMS.  INC.  SAN  DIEGO.  C A (SNODGRASS) 

(XTEAN  ENGINEERS SAUSALITO.CA(RYNECKI) 

(XTEAN  RESOURCE  ENG.  INC.  HOUSTON. TX(ANDERSON) 

OFFSHORE  DEVELOPMENT  ENG.  INC.  BERKELEY. CA 
PACIFIC  MARINE  TECHNOLOGY  LONG  BEACH.  CA  (WAGNER) 

PORTLAND  CEMENT  ASSOC.  Skokie  IL  (R-sch  & Dev  Lab.  Lib.) 

RAYMOND  INTERN  ATIONAL  INC.  E Colie  Soil  Tech  Dept.  Pennsauken.  NJ 
SAFETY  SERVICES.  INC.  A.  Patton.  Providence  RI 

SANDIA  LABORATORIES  Albuquerque.  NM  (Vortman);  Library  Div..  Livermore  CA 
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SCHUPACK  ASSOC  so.  NORWALK.  CT(SCHUPACK) 

SEATECH  CORP.  MIAMI.  FL  (PERONI) 

SHELL  DEVELOPMENT  CO.  HOUSTON..  TX  (TELES):  Houston  TX  (C.  Sellars  Jr.) 

SHELL  OIL  CO.  Houston  TX  (R.  de  Castongrene) 

SWEDEN  VBB (Library).  Stockholm 

TEXTRON  INC  BUFFALO.  NY  (RESEARCH  CENTER  LIB.) 

TIDEWATER  CONSTR.  CO  Norfolk  VA  (Fowler) 

TRW  SYSTEMS  CLEVELAND.  OH  (ENG.  LIB.);  REDONDO  BEACH.  CA  (DAI) 

UNITED  KINGDOM  British  Embassy  (Info.  Offr).  Washington  DC:  Cement  & Concrete  Assoc  (G.  Somerville) 
Wexham  Springs.  Slou;  D.  New.  G.  Maunsell  & Partners.  London;  J.  Derrington.  London;  Library.  Bristol:  R. 
Rudham  Oxfordshire:  Taylor.  Woodrow  Constr  (014P).  Southall.  Middlesex;  Taylor.  Woodrow  Constr (Stubbs). 
Southall.  Middlesex;  Univ.  of  Bristol  (R.  Morgan).  Bristol 
UNITED  TECHNOLOGIES  Windsor  Locks  CT(  Hamilton  Std  Div..  Library) 

WATT  BRIAN  ASSOC  INC.  Houston.  TX 

WESTINGHOUSE  ELECTRIC  CORP.  Annapolis  MD(Oceanic  Div  Lib.  Bryan);  Library.  Pittsburgh  PA 
WEYERH  AEUSER  CO.  LONGVIEW.  WA  (TECH  CTR  LIB) 

WISS.  JANNEY.  ELSTNER.  & ASSOC  Northbrook.  IL(J.  Hanson) 

WM  CLAPP  LABS-BATTELLE  DUXBURY.  MA  (LIBRARY) 

WOODWARD-CLYDE  CONSULTANTS  PLYMOUTH  MEETING  PA  (CROSS.  Ill) 

ADAMS.  CAPT(RET)  Irvine.  CA 
BRAHTZ  La  Jolla.  CA 

BRYANT  ROSE  Johnson  Div.  UOP.  Glendora  CA 

BULLOCK  La  Canada 

KRUZIC.T.P.  Silver  Spring.  MD 

CAPT  M U RPH  Y Sunnyvale.  C A 

GREG  PAGE  EUGENE.  OR 

R.F.  BESIER  Old  Saybrook  CT 

R.Q.  PALMER  Kaitua.  HI 

T.W.  MERM  EL  Washington  DC 

CEC  Donofrio.  John  L..  LT 


